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ABSTRACT 

The  U.S.  Marine  Corps  deploys  Marine  Air-Ground  Task 
Forces  (MAGTFs)  by  airlift  or  sealift  to  participate  in 
numerous  short-term  exercises.  These  exercises  are  of  such 
duration  that  resupply  of  the  MAGTF  by  strategic  airlift  or 
sealift  is  not  practical.  Thus,  only  stocked  spare  parts 
are  available  for  repairs  during  the  exercise. 

A  model  is  developed  which  provides  the  operational 
commander  with  a  stockage  policy  for  spare  secondary 
reparables  (e.g.,  tank  engines,  amtrack  transmissions,  etc.) 
that  optimizes  the  probability  of  successful  mission 
completion  subject  to  weight  or  volume  constraints  imposed 
by  the  MAGTF 's  mode  of  deployment.  Optimization  of  this 
stockage  policy  is  stochastically  modeled  using  data  from 
the  Marine  Corps  Integrated  Maintenance  Management  System 
data  base  and  then  solved  as  an  integer  program. 

The  integer  program  is  coded  using  the  Generalized 
Algebraic  Modeling  System  language  and  solved  using  the 
Zero/One  Optimization  Methods  mixed  integer  program  solver. 
Operational  data  for  a  Marine  Amphibious  Unit  yields  an 
integer  program  with  190  binary  variables  and  26 
constraints.  A  solution  within  0.07%  of  optimality  is 
obtained  on  an  IBM  3033AP  computer  in  3.9  seconds  and  on  a 
Zenith  Z-248  personal  computer  in  176  seconds. 
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I.   INTRODUCTION 

The  combat  structure  of  the  United  States  Marine  Corps 
(USMC)  centers  around  the  concept  of  the  Marine  Air-Ground 
Task  Force  (MAGTF)  which  is  a  combined  arms  force  consisting 
of  ground,  air,  and  combat  service  support  forces.  The  size 
of  a  MAGTF  and  its  specific  composition  are  dependent  upon 
the  particular  mission  assigned  and  the  capability  of  the 
opposing  forces.  However,  all  MAGTFs  are  organized  to  take 
maximum  advantage  of  the  combat  potential  inherent  in  a 
rapidly  deployable,  closely  integrated  air-ground  team  under 
the  control  of  a  single  commander.  A  MAGTF  consists  of  four 
basis  elements:  command  element,  ground  combat  element, 
aviation  combat  element  and  combat  service  support  element. 
It  is  the  mission  of  the  combat  service  support  element  to 
provide  all  logistical  support  for  the  ground  and  aviation 
combat  elements.  [Ref.  l:pp. 1-4]. 

In  this  study,  a  model  is  developed  which  provides  the 
operational  commander  with  a  stockage  policy  for  spare 
secondary  reparables  (e.g.,  tank  engines,  amtrack 
transmissions,  etc.)  that  optimizes  the  probability  of 
successful  mission  completion  within  the  constraints  imposed 
by  the  MAGTF' s  mode  of  deployment  (i.e.,  airlift  or 
sealift).   The  problem  of  determining  this  optimal  stockage 


policy  is  stochastically  modeled  using  operational  data 
resident  in  the  Marine  Corps  Integrated  Maintenance 
Management  System  (MIMMS)  data  base  and  then  solved  as  a 
mathematical  program.  The  remainder  of  this  chapter  will 
include  a  description  of  the  scenarios  which  are  of 
particular  interest  to  the  U.S.  Marine  Corps,  an  explanation 
of  terms,  and  a  definition  of  the  problem.  Additionally, 
^  justification    for   the    approach   taken   here    as    compared    to 

r 

L  those  in  standard  inventory  models  is  provided. 

2  A.   PROBLEM  STATEMENT 

Lu  The  U.S.  Marine  Corps  currently  deploys  MAGTFs  by 

0 

o  airlift  or  sealift  to  participate  in  numerous  short-term 

exercises  such  as  BRIGHT  STAR  in  Egypt,  GALLANT  EAGLE  in  the 
southwestern  United  States,  and  COLD  WINTER  in  Norway. 
These  MAGTFs  range  in  size  from  2,000  man  Marine  Amphibious 
Units  ( MAUs )  to  16,000  man  Marine  Amphibious  Brigades 
(MABs),  depending  upon  the  particular  exercise.  These 
exercises  are  of  such  duration  that  resupply  of  the  MAGTF's 
ground  combat  element  by  strategic  airlift  or  sealift  from 
logistics  bases  in  the  continental  United  States  is  not 
practical.  Furthermore,  by  doctrine  [Ref.  l:pp.  1-4],  a 
MAGTF  is  task  organized  to  maximize  combat  power;  therefore, 
the  number  of  mechanics,  technicians,  and  other  logistics 
personnel  together  with  their  test  and  repair  equipment  are 
kept  to  a  minimum.    Thus,  the  quantities  of  spare  repair 


parts  stocked  by  the  combat  service  support  elements  are  all 
that  are  available  to  support  the  requirements  of  the  ground 
combat  element  for  the  duration  of  the  exercise. 

In  addition  to  these  operational  exercises,  the  Marine 
Corps  deploys  Marine  Amphibious  Units  (MAUs)  to  the 
Mediterranean  Sea  as  part  of  the  U.S.  Sixth  Fleet  on  a 
continuous,  six  month  rotating  basis.  These  Mediterranean 
MAUs  remain  aboard  ship  except  when  conducting  exercises. 
Therefore,  the  inventory  of  spares  maintained  by  their 
combat  service  support  element,  called  a  MAU  Service  Support 
Group  (MSSG),  is  limited  by  the  cargo  capacity  of  the 
assigned  amphibious  ships.  Even  though  resupply  from  the 
continental  United  States  by  air  is  possible,  it  is  desired 
to  reduce  dependency  on  this  costly  method  [Ref.  2: p.  1], 
Thus,  similar  to  the  above  exercise  scenarios,  the  MSSG  must 
stock  sufficient  spares,  within  the  capacity  constraints  of 
the  assigned  shipping,  so  as  to  satisfy  the  anticipated 
demands  generated  by  the  ground  combat  element's 
participation  in  exercises  during  the  course  of  the  six 
month  deployment. 

Presently,  a  problem  of  even  greater  concern  to  Marine 
Corps  combat  service  support  organizations  than  either  of 
those  previously  mentioned  is  the  Marine  Corps's  role  in  the 
rapidly  evolving  maritime  strategy  of  the  U.S.  Navy.  As  the 
amphibious  power  projection  arm  of  the  Navy,  the  Marine 


Corps,  in  conjunction  with  the  Navy,  has  developed  and 
recently  implemented  the  concept  known  as  the  Maritime 
Prepositioning  Force  (MPF).  According  to  doctrine  [Ref.  3], 
a  MPF  consists  of  a  squadron  of  four  to  five  ships  loaded 
with  combat  equipment  and  supplies.  One  such  MPF  squadron 
is  presently  located  in  each  of  the  following  three  ocean 
areas:  Eastern  Atlantic,  Indian  Ocean, and  Western  Pacific. 
When  a  MPF  operation  is  ordered,  a  fly-in  echelon  consisting 
of  16,500  Marines  and  sailors,  comprising  a  Marine 
Amphibious  Brigade  ( MAB )  and  a  Navy  Support  Element,  is 
airlifted  to  a  benign  airfield  in  the  vicinity  of  the 
objective  area  for  link-up  with  equipment  and  supplies 
aboard  the  MPF  squadron  ships.  The  Brigade  can  be  combat- 
capable  and  ready  to  move  to  an  objective  in  five  days  or 
less  and  can  operate  for  thirty  days,  independent  of  any 
strategic  airlift  or  sealift  from  the  United  States,  using 
the  supplies  and  spare  parts  inventories  on  the  ships  [Ref. 
4:p.  11].  It  is  this  last  requirement,  to  operate 
independently  for  thirty  days  using  organic  supplies  and 
inventories,  that  poses  the  greatest  challenge  to  Marine 
Corps  combat  service  support  planners.  However,  as  with 
supporting  short-term  exercises  and  MAU  deployments,  the 
problem  of  supporting  the  MPF  MAB  simplifies  to  one  of 
stocking  adequate  spare  parts  inventories,  subject  to  the 
capacity  constraints  of  the  MPF  squadron  ships,  to  meet  the 


demand  during  the  MAB ' s  projected  missions.  The  fact  that 
such  demand  may  vary  unpredictably  must  be  confronted. 

At  the  present  time,  to  solve  spare  stockage  problems 
similar  to  those  described  above.  Marine  Corps  combat 
service  support  planners  determine  the  number  of  secondary 
items  (referred  to  in  Marine  Corps'  terminology  as  secondary 
reparahles )  to  stock  for  each  principal  end  item  in  the 
particular  MAGTF's  Table  of  Equipment  (T/E)  by  using  either 
last-period  demand  method  or  a  simple  moving  average  of  past 
or  historical  demands  generated  during  similar  exercises  or 
Mediterranean  deployments.  In  the  situation  where  no 
historical  demands  are  available,  such  as  with  the  new  MPF 
MABs,  stockage  levels  are  determined  using  estimated 
replacement  rates  in  combination  with  the  experience  of 
maintenance  and  supply  personnel.  These  techniques  have 
proved  inadequate  [Ref.  2:p.  1]  because  they  fail  to 
consider  such  factors  as  changes  to  a  MAGTF's  Table  of 
Equipment  and  variations  in  the  operational  schedule,  both 
of  which  may  occur  from  one  exercise  or  deployment  to  the 
next.  More  importantly,  they  neglect  to  consider  that  the 
demands  placed  on  the  supply  system  are  the  direct  result  of 
equipment  failures. 

To  overcome  these  inadequacies,  it  is  the  objective  of 
this  study  to  develop  a  model  which  provides  the  operational 
commander  at  the  MAU  or  MAB  level  with  a  stockage  policy  for 


spare  secondary  reparables  that  optimizes  the  probability  of 
successful  mission  completion  within  the  constraints  imposed 
by  the  mode  of  deployment.  The  proposed  model  is 
sufficiently  general  so  as  to  be  adaptable  to  any  of  the 
above  scenarios  and  it  can  accommodate  any  additions  or 
deletions  to  the  MAGTF's  Table  of  Equipment  or  variations 
in  mission  duration  among  any  combination  of  principal  end 
items.  Additionally,  the  model  incorporates  a  measure  of 
effectiveness  which  provides  the  operational  commander  with 
an  indicator  of  the  overall  impact  of  the  stockage  policy  on 
successful  mission  completion.    Furthermore,  the  model 
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Marine  Corps  Integrated  Maintenance  Management  System 
(MIMMS)  to  determine  operationally  based  failure  rates. 
Thus,    by   use    of    this    model    a    nearly    optimal    stockage    policy 
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determined  for  a  wide  variety  of  inventory  problems 
currently  facing  Marine  Corps  combat  service  support 
planners.  This  will  result  in  a  more  effective  utilization 
of  the  limited  space  and  weight  resources  available  for 
spares'  inventories  on  shipping  or  aircraft  assigned  to 
deploying  Marine  Corps  units. 


B.   INVENTORY  MODELS  AND  SOLUTION  APPROACHES 

The  stockage  problems  confronting  Marine  Corps  planners 
can  all  be  viewed  as  variations  of  the  single  period, 
constrained,  multi-item  inventory  problem.  The  simplest 
such  single-period  model  is  the  classic  newsboy  problem 
[Ref.  5],  alternatively  referred  to  as  the  Christmas  tree 
problem.  The  characteristic  distinguishing  the  models  of 
this  thesis  from  other  standard  inventory  models,  such  as 
the  classical  economic-order-quantity  or  economic-order- 
interval  models  [Ref.  6:pp.  159-289]  is  that  only  a  single 
time  period,  usually  of  finite  length,  is  relevant; 
therefore,  only  a  single  initial  procurement  or  order  is 
made  [Ref.  5:p.  297].  The  previously  described  scenarios 
all  fall  within  the  context  of  single-period  models  because 
the  demand  for  secondary  reparables  is  the  direct  result  of 
equipment  failures  which  occur  at  infrequent  intervals 
during  the  course  of  the  exercise.  Resupply  either  does  not 
occur  or  is  kept  to  a  minimum,  and  thus  it  can  be  ignored. 
Furthermore,  it  is  Marine  Corps  policy  to  permit  only 
functional  unit  replacement  at  the  field  or  organizational 
levels  of  maintenance.  Actual  repair  of  secondary  reparables 
at  the  intermediate  maintenance  level  is  very  limited  and 
for  all  practical  purposes  performed  only  at  the  depot 
level.  Therefore,  since  the  MAGTFs  of  concern  in  these 
scenarios  are  usually  only  authorized  up  to  the  intermediate 


level  of  maintenance,  the  actual  repair  and  return  of 
secondary  reparables  to  these  retail  levels  of  inventory  is 
insignificant  and  can  be  ignored.  This  is  particularly  true 
when  one  considers  that  some  of  the  exercise  durations  are 
of  the  same  order  of  magnitude  as  the  average  repair  time 
for  some  of  the  secondary  reparables.  Thus,  each  scenario 
consists  of  an  exercise  (or  exercises)  which  is  (are)  of 
finite  duration  and  during  which  repair  and  resupply  are 
essentially  nonexistent,  so  the  combat  service  support 
planner  must  make  a  one-time  decision  of  what  to  stock,  not 
how  much  of  each  item  to  procure.  He  must  do  so  within  the 
constraints  of  the  mode  of  deployment.  This  is  analogous  to 
the  single  period,  constrained,  multi-item  inventory 
problem. 

There  are  several  examples  in  the  literature  of  single 
period  inventory  models  applied  to  military  logistics 
situations.  They  include  the  spare  parts  kit  problem  [Ref. 
7:pp.  281-295],  the  fly-away  kit  problem  [Ref.  6:p.  324], 
and  the  submarine  provisioning  problem  [Refs.  8: pp.  235-243, 
9:p.  407],  The  problem  with  each  of  these  military 
applications  is  that  they  assume  that  there  is  some  monetary 
stockout  cost  incurred  when  a  demand  exists  and  no  spares 
are  available  to  meet  the  demand.  This  is  realistic  in  the 
civilian  sector  where  real  monetary  values  can  often  be 
associated  with  each  stockout.   But,  in   most  practical 


8 


military  applications  it  is  difficult,  if  not  impossible,  to 
determine  a  tangible,  realistic,  and  objective  stockout 
cost.  This  is  particularly  true  for  each  of  the  previously 
described  scenarios  because  there  is  no  way  to  place  a  cost 
on  the  inability  to  complete  a  mission  because  of  a 
deficiency  of  spares.  Thus,  it  is  difficult  to  apply  the 
models  existing  in  the  literature  to  these  real-world  Marine 
Corps  inventory  problems. 

This  problem  is  not  insurmountable.  In  the  civilian 
sector  when  an  organization  does  not  know  its  stockout  costs 
or  is  not  confident  in  its  estimates  of  them,  management 
will  set  service  levels  which  indicate  the  ability  to  meet 
customer  demands  from  stock  [Ref .  10: p.  149] .  A  service 
level  is  basically  the  complement  of  the  probability  of 
stockout  and  really  represents  a  subjective  decision  by 
management  to  accept  some  degree  of  stockout  or  customer 
disservice  [Ref.  10:p.  150].  The  U.S.  Navy  actually  uses  a 
similar  concept  when  establishing  its  Aviation  Consolidated 
Allowance  Lists  (AVCAL)  and  Coordinated  Shipboard  Allowance 
Lists  (COSAL)  [Ref.  ll:pp.  4-30-4-50].  The  service  level 
measure  is  less  than  satisfactory  for  military  applications 
since  it  is  subjectively  established  by  combat  service 
support  planners.  Service  levels  do  not  provide  the 
operational  commander  in  the  field  with  a  meaningful 
indicator  of  the  impact  of  the  recommended  stockage  policy 
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on  successful  mission  completion.  This  is  because  there  is 
no  direct  relationship  between  the  probability  of  stockout 
and  the  probability  of  successfully  completing  a  given 
mission.  A  high  service  level  does  not  equate  to  an  optimal 
probability  of  mission  completion.  This  is  particularly 
true  in  the  constrained,  multi-item  inventory  situation 
where  it  might  actually  be  better  to  accept  a  greater 
probability  of  stockout  (lower  service  level)  for  a 
secondary  reparable  with  a  low  failure  rate  in  order  to  be 
able  to  stock  more  of  another  secondary  reparable  with  a 
high  failure  rate. 

To  overcome  the  disadvantages  associated  with  the 
^  minimization  of  stockout  costs  (actually  variable  costs)  and 

ui  service  level  approaches,  the  proDosed  model  introduces  a 

u 

D 

Q  new  measure  of  effectiveness,  the  probability   of  mission 

0  —  

a: 

<^  completion     given    the    specified    stockage   policy ,    which  will 

be  referred  to  as  the  probability    of   mission    completion 

throughout  the  remainder  of  this  paper.    This  probability 

extends  the  familiar  concept  of  reliability  to  include  a 

system  and  its  inventory  of  spare  parts  [Ref.  8: p.  237].   A 

system  is  defined  to  be  the  quantity  of  a  principal  end  item 

required  to  be  operational  throughout  the  duration  of  a 

given  mission.    This  end  item  must  be  a   "readiness 

reportable  pacing  item"  as  delineated  in  Marine  Corps 

Bulletin  3000  [Ref.  12:pp.  1-12  End  (2)].   A  readiness 
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reportable  pacing  item  is  a  principal  end  item  whose 

operation  is  essential  to  the  ability  of  a  MAGTF  to  complete 

its  assigned  mission.  Simply  stated,  the  probability  of 

mission  completion  is  the  probability  that  k-out-of-n  (k£n) 

readiness  reportable  pacing  items  of  the  same  type  remain 

operational  for  a  specified  mission  with  a  given  stockage 

policy.   Black  and  Proschan  [Ref.  7:p.  283]  proposed  a 

similar  measure  of  effectiveness  in  their  early  work  using 

marginal  analysis  to  determine  optimal  spare  parts  kits  at 

minimum  cost.    They  called  it  probability    of    adequacy    or 

assurance  of  adequacy ,    which  is  defined  as 

The  probability  that  either  the  equipment  does  not  fail 
due  to  a.  .  .malfunction  or,  if  it  does  fail  for  this 
reason,  a  replacement  is  available  to  provide  continued 
successful  equipment  operation.   [Ref.  8:p.  237] 

However,  this  definition  applies  only  to  series  systems  in 

which  all  equipment  must  remain  operational  for  a  successful 

mission.   The  proposed  probability  of  mission  completion  is 

more  general  since  it  can  be  extended  to  "k-out-of-n 

systems".    This   is  required  in  the  present  study  because 

the  operational  commander,   not  the  military  logistician, 

determines  the  quantities  of  a  principal  end  item  which  he 

feels  are  necessary  to  complete  the  assigned  mission.   For 

example,  a  MAU  commander  may  feel  that  he  can  accomplish  a 

particular  mission  with  just  three  out  of  five  tanks  which 

permits  fewer  tank  secondary  reparables  to  be  stocked  than 

if  he  required  all  five  tanks  to  complete  the  same  mission. 

11 


In  this  study  the  inventory  problems  in  each  of  the 
previously  described  scenarios  are  modeled  as  single  period, 
constrained,  multi-item  inventory  problems  in  which  the 
demands  are  the  result  of  independent  failures  of  secondary 
reparables.  In  the  mathematical  formulation  of  these 
models,  the  probability  of  mission  completion  is  maximized, 
as  opposed  to  the  standard  approach  of  minimizing  expected 
(variable)  costs.  This  provides  models  which  reflect  the 
fact  that  in  most  military  situations,  and  specifically  in 
the  three  Marine  Corps  scenarios  studied  here,  stockout  or 
shortage  costs  cannot  be  determined.  The  inclusion  and 
numbers  of  constraints  are  dependent  on  the  mode  of 
deployment.  However,  in  all  cases  studied,  they  are  linear 
and  will  be  referred  to  as  capacity  constraints.  Thus,  each 
model  will  be  viewed  as  a  mathematical  program  in  which  the 
objective  function  to  be  maximized  is  a  nonlinear, 
nondecreasing  function  of  the  decision  variable  representing 
the  quantity  of  secondary  reparables  to  stock.  Inclusion  of 
the  linear  capacity  constraint  or  constraints  result  in  an 
overall  integer  nonlinear  programming  problem.  Since  the 
objective  function  is  itself  a  complicated  function  of  the 
decision  variable,  the  solution  approach  is  to  linearize  the 
objective  function  and  then  reformulate  the  entire  problem 
using  binary  variables.  This  reformulation  converts  the 
original  integer  nonlinear  programming  problem  into  a  binary 
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linear  programming  problem  which  is  then  coded  for  computer 
processing  using  the  Generalized  Algebraic  Modeling  System 
(GAMS)  [Refs.  13,  14,  15]  and  solved  using  the  Zero/One 
Optimization  Methods  (ZOOM)  [Ref.  16]  mixed  integer  program 
solver. 

Since  one  of  the  primary  purposes  of  this  study  is  to 
provide  the  combat  service  support  planners  of  the  Second 
Force  Service  Support  Group  at  Camp  Lejeune,  North  Carolina 
with  models  solvable  on  a  personal  computer,  it  is  desired 
to  avoid  dependency  on  commercial  optimization  packages  not 
locally  available.  Therefore,  the  initial  solution 
methodology  utilized  Langrangian  relaxation  to  establish  an 
upper  bound  on  the  optimal  solution  of  the  integer  linear 
problem.  Then  a  heuristic  based  on  marginal  analysis  is 
used  to  improve  upon  the  best  feasible  solution  obtained 
during  the  establishment  of  this  upper  bound,  and  this 
heuristically  improved  solution  is  accepted  as  the  final 
solution  provided  that  its  associated  objective  value  was 
reasonably  close  to  the  upper  bound.  This  methodology  is 
similar  to  that  implemented  by  DeWolfe  [Ref.  17:pp.  28-32] 
to  solve  a  Selective  Reenlistment  Bonus  problem  for  the  U.S. 
Marine  Corps.  The  advantage  of  this  approach  is  that  it  can 
be  programmed  in  a  language  such  as  FORTRAN  77  thus  avoiding 
dependency  on  a  commercial  optimizer.  However,  as  the  number 
of   constraints   increase   the   programming   difficulty 
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associated  with  the  search  procedure  used  in  establishing 
the  best  upper  bound  increases,  and  for  any  more  than  two 
constraints  actually  becomes  prohibitive.  So,  to  keep  the 
solution  approach  as  flexible  and  general  as  possible,  it 
was  decided  to  solve  all  integer  linear  problems  with  the 
GAMS/ZOOM  optimization  package  ( commerically  available  in 
March  1988),  thus  allowing  the  solution  of  models  with  more 
than  one  constraint  and  avoiding  restrictive  programming 
considerations . 

The  above  solution  approach  overcomes  the  disadvantages 
associated  with  more  standard  solutions  of  the  single  period 
inventory  model.  For  instance,  it  could  be  argued  that 
since  these  inventory  problems  are  all  examples  of  decision 
making  under  risk,  they  could  be  solved  using  decision 
theory  and  payoff  matrices  to  minimize  expected  costs. 
However,  as  was  pointed  out  earlier,  realistic  monetary 
stockout  costs  are  extremely  difficult  to  determine  in  most 
military  applications.  Additionally,  in  the  multiple  item 
inventory  problem,  the  payoff  matrices  become  very  large  and 
the  interrelationship  among  the  payoffs  associated  with 
different  combinations  of  items  becomes  difficult  to 
evaluate  and  enumerate,  so  it  is  computationally  more 
efficient  to  use  a  marginal  analysis  approach  [Ref.  10:p. 
269].  Marginal  analysis  provides  the  basis  for  the 
techniques  used  in  the  spare  parts  kit  problem  [Ref. 
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7 : pp. 283-284]  and  in  the  submarine  provisioning  problem 
[Ref.  8:p.  237],  but  it  has  practical  limitations  when  more 
than  one  constraint  must  be  considered.  Hadley  and  Whitin 
[Ref.  6 : pp  328-330]  obtained  exact  integer  solutions  to 
their  single  constraint,  fly-away  kit  problem  using  dynamic 
programming  and  Dreyfus  and  Law  [Ref.  18:pp.  107-118]  offer 
dynamic  programming  algorithms  to  solve  single  constraint, 
cargo-loading  problems.  However,  they  all  point  out  that 
for  problems  involving  more  than  two  constraints  even 
dynamic  programming  has  computational  limitations, 
particularly  when  both  the  number  of  constraints  and 
decision  variables  are  large  [Refs.  6:p.  331,  18:pp.  107- 
118]  .  The  solution  approach  utilized  in  this  study  avoids 
these  shortcomings. 

C.   THESIS  OUTLINE 

This  thesis  develops  and  presents  a  method  for 
determining  spare  stockage  levels  of  secondary  reparables 
associated  with  the  readiness  reportable  principal  end  items 
of  a  given  Marine  Air-Ground  Task  Force.  In  Chapter  II,  the 
maintenance  data  base,  data  analysis,  and  problems 
encountered  during  this  analysis  are  discussed.  The 
objective  functions  are  defined,  developed,  and  justified  in 
Chapter  III.  Additionally,  each  scenario  is  modeled  and 
formulated  as  an  integer  nonlinear  programming  problem  for 
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which  the  objective  function  is  linearized  so  that  the 
problem  can  be  reformulated  as  a  integer  linear  problem  to 
facilitate  its  solution.  Specific  solution  methodology  is 
presented  in  Chapter  IV.  The  initial  approach  using 
Lagrangian  relaxation  with  a  heuristic  is  briefly  described 
before  the  more  general  procedure  using  GAMS/ZOOM  is 
developed  and  implemented.  In  Chapter  V  results  using 
actual  data  for  a  typical  Marine  Amphibious  Unit  are 
u  presented  and  verification  of  the  model  is  conducted  using 

f  different  sets  of  simulated  failure  data.  The  chapter 

^  concludes  with  a  discussion  of  computational  experience  for 

ui 

Q  both  the  mainframe  and  personal  computer.    Chapter  VI 

I-  contains  conclusions  and  recommendations.   Finally,  listings 

2  of  both  FORTRAN  77  and  GAMS  source  codes  are  provided  in  the 

u 

Q  appendices. 
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II,   DATA  ISSUES  AND  ANALYSIS 

A  primary  impetus  for  this  study  was  the  desire  to 
develop  a  reliability-based,  spares  stockage  model  for 
secondary  reparables  utilizing  the  maintenance  data  resident 
in  the  Marine  Corps  Integrated  Maintenance  Management  System 
(MIMMS).  Currently,  spare  stockage  models  within  the  Marine 
Corps  are  limited  to  those  based  on  forecasts  from  past 
demands  generated  against  the  supply  system.  This 
forecasting  approach  neglects  to  consider  that  the  source  of 
the  demands  is  the  failure  of  equipment  such  as  secondary 
reparables,  and  it  is  precisely  this  type  of  failure 
information  that  can  be  extracted  from  relevant  MIMMS  sub- 
files which  permits  the  development  of  a  reliability-based 
model. 
A.   DESCRIPTION  OF  DATA 

MIMMS  is  an  automated  system  which  accumulates 
maintenance  data  on  all  serialized,  principal  end  items  in 
the  Marine  Corps  inventory.  It  records  and  accumulates,  in 
various  sub-files,  maintenance  actions  performed  on  a 
serialized  end  item  during  its  service  life  at  each  of  the 
three  levels  of  maintenance  (i.e.,  organizational, 
intermediate,  and  depot)  [Ref.  19].  In  particular,  the 
MIMMS  sub-file  of  relevance  to  this  study  is  the  Equipment 
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Repair  Order  History  File.  This  sub-file  maintains  all 
equipment  repair  orders  ( EROs )  that  were  completed  on  every 
serialized  end  item  during  the  past  eighteen  months,  and  it 
is  updated  on  a  quarterly  basis.  An  ERO  is  opened  each  time 
a  serialized  end  item  is  inducted  into  the  maintenance  cycle 
and  contains  information  such  as  status  of  the  end  item 
(i.e.,  combat  deadlined,  safety  deadlined,  mission  degraded, 
etc.),  the  meter  reading  at  the  time  of  the  current 
u.  maintenance  action  (expressed  in  either  hours,  miles,  rounds 

~  or  days  depending  on  the  end  item),  and  the  defect  code 

which  describes  the  nature  of  the  failure.   By  looking  at 
only  those  EROs  in  the  MIMMS  ERO  History  File  which  were 
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H  opened  for  critical,  combat  deadlining  repairs  (i.e., 
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category  code  M),  a  new,  local  user  file  can  be  created, 


Q 
lij 
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Q  This  user  file  is  organized  by  principal  end  item  and 

0 
o: 

c-  includes   only   those   end   items   which   are   considered 

111  -^ 

a: 

"readiness   reportable"   [Ref.   12:pp.   End   (1)   1-23]. 

Associated  with  each  end  item  are  all  serial  numbers  for 

which  a  combat  deadlining  ERO  exists;  each  serial  number  has 

a  meter  reading  and  defect  code  to  reflect  each  ERO 

completed  on  that  particular  serial  number  during  the 

eighteen  month  period  covered  by  the  ERO  History  File. 

Obviously,  there  may  be  multiple  entries  under  the  same 

serial  number  indicating  that  this  end  item  suffered  more 

than  one  combat  deadlining  failure  during  the  reporting 
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period.  There  may  be  no  entries  if  no  failures  were 
experienced  by  the  particular  serial  number.  Table  1 
provides  a  sample  extract  from  this  user-created  data  file 
and  displays  its  contents  for  several  tank  serial  numbers. 


TABLE  1 
EXTRACT  FROM  LOCAL  USER  DATA  FILE 


SERIAL  # 


EOTC* 


METER 
READING 


DEFECT 
CODE** 


502814 
502828 


502831 
502835 


H 
H 
H 
H 
H 
H 
H 
H 
H 


74 

51 

75 

107 

155 

163 

81 

149 

62 


D05  (POWER  TRAIN, CLUTCH) 

A34  (ENGINE   REPLACE) 

J34  (COOLING  SYS  REPLACE) 

C34  (POWER  RACK  REPLACE) 

A34  (ENGINE  REPLACE) 

B34  (TRANSMISSION  REPLACE) 

A34  (ENGINE  REPLACE) 

B34  (TRANSMISSION  REPLACE) 

KOI  (ELEC  SYS,  GENERATOR) 


*   EOTC=   Equipment  Operating  Time  Code 

(H=Hours,  M=Miles,  R=Rounds,  D=Days) 

**  Defect  Codes  are  alphanumeric  codes  and  have  been 
annotated  here  for  display  purposes  only. 


The  data  in  this  consolidated  file  can  be  used  to  estimate 
failure  rates  for  each  secondary  reparable  in  either  of  two 
ways  and  both  methods  will  be  described  in  the  following 
section  along  with  their  advantages  and  disadvantages. 
It  must  be  stated  that  both  the  format  and  method  of  update 
for  many  of  the  MIMMS  files  require  improvement  before  a 
truly  comprehensive  analysis  of  the  data  is  possible.  A 
case  in  point  is  the  ERO  History  File  used  in  this  study, 
and  recommendations  for  its  improvement  are  provided  in  the 
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following  paragraph  and  in  Chapter  VI  which  minimize  the 

impact  on  current  MIMMS  file  management  and  the  amount  of 

effort   required  at  the  using  unit  input  level. 

One  of  the  basic  problems  with  the  ERO  History  File  is 

that  potentially  useful  data  is  being  discarded  during  each 

quarterly  update  in  order  to  minimize  file  size.   What  is 

proposed  here  is  to  maintain  a  local   (i.e.,  at  the  Force 

Service  Support  Group  level )  consolidated  version  of  the  ERO 

History  File  with  a  format  similar  to  that  described  in  the 

preceding  paragraph.   This  file  should  be  updated  when  its 

parent  ERO  History  File  is  updated  at  the  end  of  each 

quarter.    Since  this  consolidated  file  contains  only 

information  extracted  from  combat  deadlining  EROs,  which 

comprise  a  small  fraction  of  the  total  EROs  in  the  parent 

file,  there  is  a  greater  capability  to  accumulate  relevant 

data  without  the  need  to  automatically  discard  the  oldest 

three  months  of  data.    The  basic  algorithm  for  each 

quarterly  update  is: 

Step  1.  Select  a  combat  deadlining  ERO  (Category  Code  M) 
from  the  ERO  History  File  if  one  exists.  If  none 
exist,  STOP. 

Step  2.  Determine  if  the  serial  number  associated  with 
this  ERO  is  already  resident  on  the  consolidated 
file.  If  it  is,  continue  with  Step  3.  If  it  is 
not,  add  the  serial  number  and  its  associated 
meter  reading  and  defect  code  to  the  consolidated 
file  under  the  appropriate  principal  end  item. 
Return  to  Step  1. 
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step  3.  Determine  if  the  same  defect  code  exists  for  this 
serial  number.  If  it  does,  continue  with  Step  4. 
If  it  does  not,  add  the  meter  reading  and  defect 
code  to  the  list  associated  with  this  serial 
number.  Return  to  Step  1. 

Step  4.  Subtract  the  two  meter  readings  to  obtain  a  time 
between  failure  for  this  type  of  defect.  Add 
this  new  data  point  to  the  list  associated  with 
this  type  of  defect.  Delete  the  record 
associated  with  the  older  meter  reading  (i.e., 
the  smallest  meter  reading)  and  add  the  new  meter 
reading  and  defect  code  to  the  end  of  the  list 
associated  with  this  serial  number.  Return  to 
Step  1. 

The  result  is  an  updated  version  of  the  consolidated  file 

which  is  retained  at  the  local  level  for  update  during  the 

next  quarterly  cycle.   Each  readiness  reportable  principal 

end  item  in  the  file  is  broken  down  by  serial  number  which 

in  turn  contain  records  of  the  most  recent  meter  readings 

and  defect  codes.   Each  end  item  also  consists  of  lists 

containing  times  between  failures  for  each  type  of  defect. 

Thus  this  local  consolidated  file  does  not  alter  any 

existing  MIMMS  file  definitions,  nor  does  it  necessitate  the 

creation  of   any  new  system   file   definitions.     More 

importantly,  it  provides  the  combat  service  support  planner 

with  a  consolidation  of  all  the  pertinent  information 

required  to  compute  failure  rates  and  track  the  life  history 

of  any  secondary  reparable  associated  with  a  readiness 

reportable  end  item.   This  is  a  significant  improvement  over 

the  present  method  which  is  restricted  to  using  the  most 

recent  eighteen  months  of  failure  data  resident  in  the  ERO 


21 


Q 
Ui 

u 


History  File,  and  it  is  this  restriction  that  limits  the 
accumulation  of  failure  times  for  certain  types  of  secondary 
reparables  (e.g.,  certain  radio  transmitters). 

B.   DATA  ANALYSIS 

Using  the  data  resident  in  the  consolidated  file,  the 
operating  times  between  failures  for  each  type  of  secondary 
reparable  were  computed.  Failure  rates  were  estimated  using 
two  different  approaches.  In  the  first  approach,  an 
exponential  lifetime  distribution  was  assumed  for  each 
secondary  reparable  and  the  maximum  likelihood  estimator 
(M.L.E.)  for  the  parameter  ).^   of  an  exponential  distribution 


Ui 

> 
O 
^  was  used  to  provide  an  estimate  of  the  failure  rate.  The 


second  approach  involves  a  Bayesian  analysis  of  the  times 
between  failure  to  determine  a  posterior  distribution  for 
the  failure  rate  of  a  secondary  reparable. 
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0 
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c:  In  deriving  the  M.L.E.  for  X^,     it  is  assumed  that  the 


observed  lifetimes  of  each  secondary  reparable  i  form  a 
relevant  random  sample  of  size  n^^  from  an  exponential 
distribution  with  unknown  parameter  i^.  The  likelihood 
function  must  also  account  for  the  exposure  times  during 
which  no  failures  of  type  i  occurred.  If  these  terms  are 
not  considered,  the  computed  failure  rate  will  tend  to  be 
too  large  and  excessive  stockage  levels  would  result.  For 
the  n^   observed  failure  times  and  the  observed  exposure  time 
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I'l    during  which  no  failures  occurred,   the   likelihood 
function  is: 


",■ 


_;.,(V;  ,  ,,^.) 


/(l- I  //)  =  A^-exp  '•'^'1  . ...  ./^exp  ''''^  •  exp  ''•''''  =  /"^  exp     '^Tj '^ 


Maximizing   the   likelihood   function   with   respect   to 

yields  the  following  estimator  of  the  failure  rate: 
A        ri;  ri; 


'I 


I    ~      j-i 


(2.1) 


where  ■J"/  represents     the  total  exposure  time  for  secondary 
reparable  i.  [Ref.  20:pp. 282-296] 

The  Bayesian  approach  also  assumes  that  the  distribution 
of  the  lifetimes  for  each  secondary  reparable  i  is 
exponential  with  parameter  / ^;  however,  it  further  assumes 
that  the  exact  value  of  this  parameter  is  a  realization  of  a 
random  variable.  Even  though  /^  is  not  precisely  known,  it 
has  a  prior  distribution  which  is  taken  here  to  be  a 
conjugate  Gamma  distribution  with  shape  parameter 
a  ^('j.^>0)    and  scale  parameter  j]  ^{   /j^>0): 

r 


Mi ;  y-i ,  /?,•)  =  < 


^^'  ;:f'-^e.xp-^'^^  ;.y>o 


na,)"'    ^"^    '-^"  (2.2) 

0  ,  //  ^  0  . 

The  prior  distribution  parameters  are  estimated  using  a 
parametric  empirical  Bayes  (PEB)  approach  for  failure  rates 
as   outlined   by   Gaver    and   Lehoczky    [Ref.21:pp.    220-224].       The 
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PEB  approach  uses  the  entire  data  set  to  compute  estimates 
a  ^  and  p^  which  are  substituted  into  the  applicable 
formulas.  To  compute  the  prior  parameters  in  PEB  requires 
the  maximization  of  the  "marginal  likelihood  function" : 


where, 
m 


'0     "im-  r(ay) 


individual  source  of  data  for  component  i  (e.g., 
tank  serial  #502826) 


^(i       total  number  of  data  sources  for  component  i 
failure  rate  for  component  i  of  data  source  m 


•im 


hm      observed  exposure  time  of  component  i  on  data  source 
m 


irn 


number  of  failures  of  component  i  on  data  source  m 

^i        shape  parameter  of  Gamma  prior  distribution  of 
failure  rate  for  component  i 

Pi       scale  parameter  of  Gamma  prior  distribution  of 
failure  rate  for  component  i. 

Integration  results  in  the  following  marginal  likelihood 

function  for  secondary  reparable  i: 

which  is  the  product  of  individual  negative  binomial 
distributions  with  parameters   p  = —       and     r  =  a^- 
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The  numerical  maximization  of  this  function  subject  to  known 
upper  and  lower  bounds  on  each  a  j^  and  p^  is  discussed  in 
Chapter  IV.  The  important  result  of  this  maximization  is  an 
estimate  for  the  shape  (  ai)and  scale  (p^)  parameters 
of  the  Gamma  prior  distribution  that  incorporates  the  actual 
observed  operational  data  for  each  secondary  reparable. 

Returning  to  the  Bayesian  analysis,  the  n^  observed 
times  between  failure  and  the  observed  exposure  times  r'j 
during  which  no  failures  occurred  are  all  a  function  of 
/.J  and  give  rise  to  the  likelihood  function. 


"i 


/(,,  I ;.,)  =  ;.r^  exp-A^5'^"''v  =  ;.?exp-'-'-'  (2.5) 

where  r^  once  again  represents  the  total  exposure  time  for 
secondary  reparable  i.  We  are  now  interested  in  the 
updated  distribution  of  /.  ^  after  the  times  between  failure 
have  been  observed  so  as  to  provide  a  current  estimate  of 
future  failure  events.  This  is  the  posterior  distribution 
and  it  is  proportional  to  the  product  of  the  likelihood 
function  (2.5)  and  the  prior  distribution  (2.2): 


/(;.,  1 L-)  oc  /(4- 1 ;,)/(/,) 


Thus, 


/(//  I  li)  =   cX;.;'exp-'^-'^0(Ar^exp-^'^0 
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where  c^  is  a  constant  factor  which  depends  on  the  observed 
times  between  failure,  but  does  not  depend  on  ).  ^.  The  value 
of  this  constant  can  be  determined  using  normalization.  The 
resulting  posterior  distribution  is: 

A 

a  Gamma  distribution  with  shape  parameter   ^^  +  =^j   and   scale 

A 

parameter  t^  + /?,  .  [Ref.  20: pp.  257-280] 

The  posterior  distribution  could  also  have  been  derived 
directly  using  the  fact  that  the  family  of  Gamma 
distributions  serves  as  a  conjugate  family  of  prior 
distributions  for  samples  from  an  exponential  distribution. 
Thus,  using  the  theorem  for  Gamma  conjugate  priors  [Ref. 
20: pp.  271-272],  the  posterior  distribution  for  /  j_  given 
the  observed  times  between  failure,  is  a  Gamma  distribution 
with  shape  parameter  n^^  +  a^^  and  scale  parameter  'J'i  +  A  . 
This  agrees  with  the  result  derived  in  the  previous 
paragraph  and  indicates  that  the  posterior  distribution  can 
be  defined  knowing  only  the  number  of  observations  and  the 
sum  of  all  the  exposure  times.  This  is  an  important  result 
because  it  enables  the  posterior  distribution  to  be  updated 
very  easily  as  new  observations  are  obtained. 

Since  the  posterior  distribution  of  Xi        has  been  shown 
to  be  a  Gamma  distribution,  its  expected  value  is: 
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E[/.,-  I  data]  = 


^/  +  ^i 


^i  -  l^i 
This  expected  value  also  happens  to  be  a  "Bayes  estimator" 

of   /.  ^  when  a  "squared  error  loss  function"  is  used  [Ref. 

20:p.  227].    This  estimator  could  also  have  been  used  to 

provide  an  estimate  for  each    A  ^.    However,  the  primary 

purpose  of  the  Bayesian  analysis  was  not  to  provide  a  point 

estimate  of  }.^,    but  rather  to  provide  a  distribution  (i.e., 

posterior  distribution)  for  each  /.  ^    which  could  be  used  to 

make  probabilistic  predictions  about  demand  functions 

involving  /^ .   This  step  is  detailed  in  the  next  chapter. 
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III.   MODEL  FORMULATION  AND  DESCRIPTION 

In  this  chapter,  the  two  different  objective  functions 
considered  in  this  study  are  defined,  derived,  and 
justified.  Additionally,  each  of  the  three  scenarios 
described  in  the  introduction  is  formulated  as  an  integer 
nonlinear  program  whose  objective  function  is  subsequently 
linearized  to  create  an  integer  linear  problem. 
A.   MODEL  DEVELOPMENT 

The  following  formulation  is  developed  to  model  the 
determination  of  secondary  reparable  inventory  levels  for 
each  of  the  three  scenarios  described  in  the  introduction. 
With  the  exception  of  the  number  of  constraints,  the 
formulation  is  the  same  for  each  scenario. 

secondary  reparable  or  component 

unit  shipping  weight  of  secondary  reparable  i 
[lbs] 

IV  weight  capacity  of  aircraft  or  ship  [lbs] 
v/         unit  volume  of  secondary  reparable  [cuft] 

V  volume  capacity  of  aircraft  or  ship  [cuft] 

"/         maximum   allowable   number   of   spares   for 
component  i 


INDICES: 

DATA: 
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FUNCTIONS: 

Pii'^i)  probability  that  the  number  of  readiness 

reportable  pacing  items  ( specified  by  the 
operational  commander)  containing  secondary 
reparable  i  will  survive  a  mission  given  x^ 
spares  are  stocked 


DECISION  VARIABLE: 

Xj  Number  of  spares  to  stock  for  secondary 

reparable  i 

FORMULATION: 


max 


]~]a(-^.)  (PI) 


0<x,<u^       /=1,...,.V 
Xi    integer 


The  objective  function  is  a  mathematical  expression  of 
the  probability  that  a  quantity  (specified  by  the 
operational  commander)  of  each  readiness  reportable  pacing 
item  in  a  MAGTF ' s  Table  of  Equipment  will  survive  a  mission 
of  specified  duration  given  a  particular  stockage  policy  for 
secondary  reparables.  This  probability  was  referred  to  in 
the  introduction  as  the  probability  of  mission  completion. 


29 


In  terms  of  this  definition,  successful  mission  completion 
hinges  on  the  survival  of  the  MAGTF's  pacing  items.  Each 
pacing  item  is  composed  of  secondary  reparables.  The 
failure  of  any  one  secondary  reparable  will  cause  the  pacing 
item  to  be  combat  deadlined  and  unable  to  complete  its 
mission.  Thus,  each  pacing  item  is  viewed  as  a  system  of 
secondary  reparables  in  series.  However,  in  a  MAGTF's  Table 
of  Equipment,  there  is  not  just  one  of  each  type  of  pacing 
item,  but  rather  a  specific  quantity  of  each  type  necessary 
to  accomplish  the  assigned  mission.  This  minimum  quantity 
required  will  not  only  vary  from  one  mission  to  another,  but 
also  among  different  MAGTF  commanders  when  confronted  with 
identical  missions.  Therefore,  the  different  pacing  items 
taken  as  a  whole  do  not  constitute  a  simple  series  system  in 
which  all  components  must  survive  the  given  mission. 

In  defining  the  objective  function,  two  different 
functional  forms  were  derived  and  studied.  Both  rely  on  the 
assumptions  that  secondary  reparables  fail  at  constant  rates 
and  that  the  failures  of  individual  secondary  reparables 
are  independent  of  one  another.  The  independence  assumption 
is  certainly  true  for  the  secondary  reparables  associated 
with  different  types  of  pacing  items;  however,  it  may  not 
necessarily  be  true  for  the  secondary  reparables  of  a 
specific  pacing  item.  For  example,  the  failure  of  M60A1 
tank  engines  is  independent  of  the  failure  of  LVTP7A1 
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amtrack  engines,  but  the  failure  of  amtrack  final  drives  may 
not  be  independent  of  amtrack  transmission  failures.  These 
assumptions  together  with  the  fact  that  secondary  reparables 
from  one  type  of  pacing  item  are  not  interchangeable  with 
those  of  another  (e.g.,  a  tank  engine  cannot  be  substituted 
for  an  amtrack  engine)  make  it  possible  to  consider  only 
secondary  reparables  without  having  to  specify  the  pacing 
items  of  which  they  are  components. 

The  functional  form  of  the  first  objective  function,  in 
which  the  failure  rate  of  each  secondary  reparable  is 
assumed  to  be  known,  is  as  follows: 


m=0 


where, 

number  of  pacing  items  containing  secondary  reparable 


a, 


ki  minimum  number  of  pacing  items  containing  secondary 
reparable  i  required  to  be  operational  throughout 
given  mission  (A^  <  <i,) 

A 

//   M.L.E.  of  failure  rate  for  secondary  reparable  i 

h  mission  duration  during  which  secondary  reparable  i 
must  operate. 

This  function  is  a  direct  consequence  of  the  constant 

failure  rate  assumption  for  each  secondary  reparable. 

Underlying  this  constant  failure  rate  is  an  inherent 

exponential  lifetime  distribution  which  implies  that  the 
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number  of  failures  defines  a  Poisson  (arrival)  process  with 
exponential  interarrival  times.  The  number  of  failures  of 
secondary  reparable  i  to  occur  in  any  fixed  interval  of 
time,  such  as  a  mission  of  duration  t^,  will  have  a  Poisson 
distribution  with  mean  ^- ±^±'  However,  since  there  are  a 
total  of  a^  pacing  items  containing  secondary  reparable  i 
and  operating  during  this  mission,  the  total  number  of 
failures  will  be  the  sum  of  a^  independent,  Poisson  random 
variables  (each  with  a  mean  ±tj_)  which  is  itself  a  Poisson 
random  variable  with  a  mean  ^</ih  •  This  is  really  only 
an  approximation  since  not  all  a^  pacing  items  may  be 
required  during  the  mission.  In  which  case  the  overall 
failure  rate  would  be  smaller  implying  that  we  are  probably 
overestimating  the  number  of  spares  required.  On  the  other 
hand,  it  can  be  argued  that  with  fewer  than  a_^  items  the 
overall  failure  rate  of  the  operational  items  may  increase 
since  the  same  mission  must  be  accomplished  with  fewer 
assets.  In  any  case,  this  can  be  viewed  as  a  mild,  worst 
case  approximation. 

Therefore,  equation  (3.1)  represents  the  probability 
that  there  are  less  than  or  equal  to  some  specified  number 
of  failures  of  secondary  reparable  i  during  a  mission  of 
duration  t^.  This  specified  number  is  the  sum  of  the  number 
of  spares,  x^,  and  the  difference  between  the  number  of 
pacing  items  containing  secondary  reparable  i  (a_^j  and  the 
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required  minimum  number  of  pacing  items  containing  secondary 

reparable  i  (k.^)  .        As  an  example,  consider  the  situation 

where  the  Table  of  Equipment  of  some  MAGTF  consists  of  five 

tanks  each  containing  one  engine  (a^=5)  and  assume  that  the 

MAGTF  commander  decides  that  to  accomplish  the  assigned 

mission  he  must  have  a  minimum  of  four  tanks  operational  at 

all  times  (Jc^  =  4).   Then,  assuming  one  spare  tank  engine 

(x_£  =  l)  is  stocked,  there  can  be  no  more  than  two  engine 

failures  ( a^+x^-k^=5+l-4=2 )  to  adhere  to  the  commander's 

guidance.   It  can  be  shown  that  this  probability  is  equal  to 

the  probability  that  the  quantity  of  each  type  of  pacing 

item   (or   equivalently   all   its   component   secondary 

reparables)  remains  above  a  specified  level  Ckj^) ,     because 

the  following  relation  must  always  hold  true: 

{number  of  failures  of  secondary  reparable  i}  + 

(number  of  secondary  reparables  i  operational  arid  in 
inventory}  = 

(total  number  of  secondary  reparables  i,  ( a_^  +  x^ ) } 
Thus, 

P( (number  of  secondary  reparables  i  operational  and 
in  inventory  at  t}>^k^)  = 

P(  (a^j^+Xy }  -  (number  of  failures  of  secondary  reparable 
i  in  €}>kj_)  = 

P( (number  of  failures   of   secondary   reparable   i  in  t} 
^a^+x^-k^). 

Therefore,  the  probability  of  mission  completion  for  the 

entire  MAGTF  is  the  product  of  the  probabilities  that  the 
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individual  secondary  reparables  exceed  a  level,  k^,  given  x^ 
spares  (or  equivalently  experience  at  most  a_£+x^-k^ 
failures).  So,  the  maximum  number  of  operational  pacing 
items  of  a  specific  type  at  any  given  moment  will  be  equal 
to  the  lowest  level  of  its  component  secondary  reparables. 
It  is  this  minimum  level  that  must  exceed  the  specified 
level  Jc^,  since  if  the  minimum  level  exceeds  k_£,  all  levels 
will  exceed  7c^.  Now,  by  the  assumption  that  the  secondary 
reparable  failures  are  independent  combined .  with  the  fact 
that  if  one  type  of  pacing  item  cannot  complete  the  mission 
the  entire  MAGTF  will  not  be  able  to  complete  the  mission, 
the  probability  that  the  minimum  level  exceeds  the  specified 
level  (i.e.,  p^(x^))  can  be  expressed  as, 

v 

pM}  'P2U2)  •  -  'Px(-Xx)  =  rip,[x^)  (3.2) 

1=1 

which  accounts  for  the  functional  form  of  the  objective 
function. 

Similar  to  the  first  objective  function,  the  second 
objective  function  also  takes  advantage  of  the  assumptions 
of  constant  failure  rates  for  each  secondary  reparable  and 
independence  among  the  failures  of  different  secondary 
reparables;  however,  it  incorporates  a  more  realistic 
approach  as  to  the  nature  of  the  failure  rates.  In  the 
previous  objective  function,  the  failure  rate  Aj^,  of  each 
secondary  reparable  was  assumed  to  be  known  (actually  an 
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estimate  computed  using  the  maximum  likelihood  estimator  of 
the  failure  rate  must  realistically  be  used).  For  this 
second  objective  function,  a  Bayesian  approach  is  taken  in 
which  the  exact  value  of  ;.  ^  is  explicitly  recognized  to  be 
unknown.  The  observed  data  (i.e.,  times  to  failure)  are 
used  to  update  the  prior  information  concerning  each  failure 
rate  to  form  a  posterior  distribution  as  discussed  in  the 
chapter  on  data  analysis.  This  posterior  distribution  is 
necessary  to  derive  the  predictive  distribution  and  the 
functional  form  of  Pj^(Xj_),  since  we  are  interested  in 
computing  the  probabilty  distribution  of  the  number  of 
failures  of  secondary  reparable  i  which  occur  during  some 
mission  of  length  tj_.  By  analogy  to  the  previous  derivation 
of  p^(Xj_),  it  is  desired  to  consider  this  number  of 

failures  because  we  want  this  quantity  to  be  less  than  some 
specified  level  (i.e.,  a^+x^-k^)  so  as  to  insure  the  entire 
MAGTF  completes  the  assigned  mission.  Assuming  that  the 
failure  rate  x  j^  is  given,  the  total  number  of  failures  of 
a_£  secondary  reparables  of  type  i  in  a  mission  of  length  tj_ 
will  have  a   Poisson   distribution   with   approximate  mean 

^i  '■  i^±' 

^(Exactly /;•  failures  in  (0,r,]  |  /•)  =  '^^'''.7'  exp"""''-'^'' 

But  /^  is  not  precisely  known  and  has  the  following 
posterior  density  (as  derived  in  the  data  analysis  chapter) 
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given  the  observed  times  to  failure  and  prior  information: 

JIX.  I  data)  =  ^'  ^  ^'^  ^        ;.?  -^  ^'- '  exp-"'^^'  +  ^'^ 
r(/7,-  +  a,) 

where, 

^/   total  exposure  time  for  secondary  reparable  i 

f^i       number   of   observations   for   secondary   reparable 
i 

A 

«/   empirical  Bayes  estimate  of  shape  parameter  for  Gamma 
(  ^-I'^^i^^i     )  pi^ior  distribution 

A 

P^       empirical  Bayes  estimate  of  scale  parameter  for  Gamma 
(  ).i\^i^^i      )    prior  distribution. 

This  can  be  rewritten  as, 

j^;.,-idata)  = — ^ • 

Now,  to  calculate  the  predictive  distribution, 

(♦OO 

P{E\3.Qi\yf.  failures  in  (0,rj)  =         i^(Exactlyy;  failures  in  (0,rj  I  /,)  J\a,  I  dataW/,- 

Jo 

which  is  upon  substitution. 


/.OO 


Evaluation  of  the  integral  yields: 

nrr       .1     rr■^          ■     m    ■^^       Ha;,- +  a, +7;)    /        r,  +  g,        ^^,  +  «:/           a,//  ^/j- 

/•(Exactly  y*  failures  in  (0,rj)  = —  (  j  (  —  \ 

f^  r(/7,-  +  a,)    a^i,  +  X,  +  )i  ^  a/,  +  r,  +  ?,•  ^ 

This  is  a  negative  binomial  distribution  with  parameters 

T, + h 

P  = ;;-    and     r  =  rt^  +  a^-  . 

a,/,-  +  T,-  +  /?,. 
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Even  though  Pj^(x^)  is  defined  as  the  probability  that 
the  stockage  of  x^  spares  for  secondary  reparable  i  allows 
the  pacing  item  of  which  it  is  a  component  to  exceed  some 
minimum  level  ( fc^ ) ,  it  has  been  shown  that  this  is 
equivalent  to  the  probability  that  the  number  of  failures  of 
secondary  reparable  i  during  mission  t j_  is  less  than  or 
equal  to  aj_+Xj^-kj_  (i.e.,  P(  number  of  failures  of  secondary 
reparable  i  in  {0,tj_J<_    a^+x^-k^)).   Therefore,  the  second 


functional  form  of  p^fx^j  is: 

a;  -t-  X:  —  kj  .  A 

'r(«,  +  5,  +  m)      /        T.-rft 


iW=      Z 


a^ti  ■  ^ 


(3.3) 


PlV^lJ  /__j  A  \  A  /  V  A 

7=0        "^!  I~("/  +  ^i)  a,ti  +  T,-  +  i^i  a,t,  +  T,-  +  /) 

where, 

flj  number  of  pacing  items  containing  secondary  reparable 
i 

k;  minimum  number  of  pacing  items  containing  secondary 
reparable  i  required  to  be  operational  throughout 
given  mission  (  k j__<  a j_ ) 

/7;  number  of  observations  ( data  points )  of  times  between 
failure  for  secondary  reparable  i 

'^i       total  observed  exposure  time  for  secondary  reparable  i 

/■/  mission  duration  during  which  secondary  reparable  i 
must  operate 

a^  empirical  Bayes  estimate  of  shape  parameter  for  Gamma 
(  /,■ ;  a;  ,  /?,•   )  prior  distribution 

A 

/?^.  empirical  Bayes  estimate  of  scale  parameter  for  Gamma 
(  // ;  ^i ,  Pi      )  prior  distribution  , 

The  reasoning  behind  the  functional  form  of  equation  (3.2) 

developed  for  the  first  p^(x^)    still  holds  true  for  the 
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newly  derived  pj_(Xj^);    therefore,  the  functional  form  of  this 
objective  function  is  still  the  product  of  the  individual 

P±i^±) 's. 

The  linear  weight  capacity  restricts  the  total  shipping 
weight  (i.e.,  secondary  reparable  plus  container)  for  all 
secondary  reparables  stocked  to  be  less  than  or  equal  to  the 
weight  capacity  of  the  ship  or  airplane.  Similarly , the 
linear  constraint  on  volume  or  space  capacity  requires  that 
all  secondary  reparables  stocked  not  exceed  the  volume 
available.  The  final  constraint  in  PI  indicates  that  there 
is  an  upper  bound  on  the  number  of  spares  for  each  secondary 
reparable  that  can  be  stocked.  This  is  a  very  real 
constraint  because  there  are  only  a  finite  number  of  each 
type  of  secondary  reparable  at  the  wholesale  inventory  level 
from  which  the  retail  inventories  used  to  support  all 
exercises,  deployments,  and  contingencies  can  be 
established.  Thus,  a  combat  service  support  planner  must 
consider  the  requirements  of  the  force  as  a  whole  so  as  not 
to  commit  too  many  secondary  reparables  to  any  one 
exercise. 

The  number  of  constraints  is  dependent  on  the  scenario. 
Any  scenario  that  involves  the  airlift  of  the  ground  combat 
element  and  its  supporting  forces  into  the  exercise  or 
operation  area,  such  as  the  short-term  exercises  described 
in  the  introduction,  will  have  only  one  constraint.   This  is 
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the  constraint  on  the  allowable  cargo  load  (ACL).   Allowable 
cargo  load  is: 

The  initial  limiting  factor  in  determining  airlift 
requirements.  .  .  .  This  is  a  weight  limiting  factor 
expressed  in  terms  of  STON  [short  tons]  which  cannot  be 
exceeded  in  planning  aircraft  loads.  The  ACL  is 
established  by  the  Air  Force  for  each  type  aircraft  for  a 
specific  mission.  [Ref.  22:pp.  4-9] 

Since  the  formulation  of  problem  PI  now  reduces  to  just  one 

constraint,  it  is  recognizable  as  a  nonlinear  knapsack 

problem  due  to  the  nonlinear  objective  function.    The 

individual  secondary  reparables  represent  the  items  or 

commodities  that  must  fit  into  a  knapsack  while  the  ACL 

constraint  corresponds  to  the  knapsack's  weight  or  volume 

capacity. 

The  other  two  scenarios  described  in  the  introduction 

both  involve  two  constraints.    These  two  constraints 

represent  the  weight  and  volume  capacities  of  the  shipping 

assigned  to  perform  the  sealift.   In  the  case  of  the  MAUs 

deploying  to  the  Mediterranean  Sea,  the  secondary  reparables 

are  normally  stored  aboard  one  of  the  ships  of  the 

amphibious  task  force  such  as  a  LHA,  LPH,  LKA,  LSD,  or  LST. 

Thus,  the  constraints  reflect  the  weight  and  volume  capacity 

of  the  ship's  cargo  hold  or  compartment  which  has  been 

designated  for  secondary  reparable  storage.    On  the  other 

hand,  in  considering  secondary  reparable  support  of  the  MPF 

MAB,  the  secondary  reparables  are  spread-loaded  among 

several  of  the  MPF  squadron  ships;  therefore,  the  weight  and 
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volume  constraints  are  actually  the  combined  weight  and 
volume  capacities  of  the  cargo  holds  of  the  specified  MPF 
squadron  ships.  Since  both  of  these  sealift  scenarios 
involve  two  contraints,  problem  PI  can  no  longer  be  regarded 
as  a  nonlinear  knapsack  problem,  but  it  is  still  an  integer 
nonlinear  programming  problem. 

Regardless  of  the  number  of  constraints  or  whether 
equation  (3.1)  or  (3.3)  is  used  as  the  functional  form  for 
p^fx^j,  the  objective  function  is  a  complex,  nonlinear 
function  of  the  decision  variable.  For  this  reason,  problem 
PI  will  be  transformed  into  a  binary  linear  programming 
problem.  The  following  section  explains  this  transformation 
and  the  new  formulation. 

B.   TRANSFORMATION  TO  A  BINARY  LINEAR  PROGRAM 

It  can  be  seen  from  the  final  constraint  of  problem  PI 
that  the  decision  variable,  Xj_,  is  constrained  to  a  set  of 
discrete  values.  Therefore  the  objective  function  will 
first  be  separated  using  the  properties  of  logarithms  and 
then  the  transformation  technique  for  converting  discrete 
variables  to  binary  variables  will  be  used  to  reformulate 
the  problem  [Ref.  23:p.  12].  Taking  the  natural  logarithm 
of  the  objective  function  yields: 
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Now,  the  final  constraint  in  PI  can  be  expressed  as, 

Xi  e  S^  =  [s^-Q,  5^1,  ...  ,  s^J  =  (0,1,  ...  ,u,} 

and  this  is  equivalent  to  the  constraint  set  [Ref.  23  :p. 
12]  : 

^  V 


S'^y  =  ^ 


/=0 


.Y 


6(0.1}    .  J  =  0,\,...,Ui 


Applying  the  separation  technique  and  transformation  to 

problem  PI  results  in  the  new  formulation: 

INDICES: 

/•  =  1 A'    secondary  reparable  or  component 

J  =  0, ...  ,u^        number  of  secondary  reparables  i 

DATA: 

fu  ln(Pi(j)) 

^ij  total  shipping  weight  for  j  spares  of  type  i 

[lbs] 

IV  weight  capacity  of  airplane  or  ship  [lbs] 

^ij  total  volume  of  j  spare  of  typei[cuft] 

V  volume  capacity  of  airplane  or  ship  [cuft] 

"/         maximum  allowable  number  of  spares  for 
component  i 

DECISION  VARIABLE: 

_  jl    if  j  spares  of  type  i  are  selected 
'j'   lO    otherwise 
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FORMULATION: 

S     "z 

max 

i=l;=0 

.V     "/ 

S.L 

E  E^^v  ^  ^^^ 

/=ly=0 

.V     "f 

E  E^V^^  ^  ^'' 

i=ly=0 

"/ 

V.V^=1            ^' 

;=0 

A>e(0,l} 

(P2) 


Assuming  known  failure  rates,  the  coefficients,  ^±-j>     i^i  the 
objective  function  are: 


fu  =  In 


where. 


;   ; exp  '  '  ^ 

i_j    m: 
L  m=0 


(3.4) 


a,-  number  of  pacing  items  containing  secondary  reparable 
i 

kj  minimum  number  of  pacing  items  containing  secondary 
reparable  i  required  to  be  operational  throughout 
given  mission  (kj^<^aj^) 

A 

>-/   M.L.E.  of  failure  rate  for  secondary  reparable  i 

^z   mission  duration  during  which  secondary  reparable  i 
must  operate. 


If  perfect  knowledge  about  the  failure  rates  is  not  assumed 
and  the  Bayesian  approach  is  adopted  the  coefficients  are: 
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a,  +;-  A:,- 


/>•  =  In 


A 


Y^  Fin.  +  7.:^  n:]      /  _JV_W^£__  \";  +  ^'  / 

L  ^0  ^'-  ^i^-i  -^  =^i)     a^ii  +  T,  4-  ^,-        ^  a.T,  +  T^-  +  /J 


m 


A 


(3.5) 


where, 

a,  number  of  pacing  items  containing  secondary  reparable 
i 

^i  minimum  number  of  pacing  items  containing  secondary 
reparable  i  required  to  be  operational  throughout 
given  mission  (  k j^_<  a^  ) 


"i 


number  of  observations  (data  points)  of  times  between 
failure  for  secondary  reparable  i 


■^/   total  observed  exposure  time  for  secondary  reparable  i 

^;    mission  duration  during  which  secondary  reparable  i 
must  operate 

QCf       empirical  Bayes  estimate  of  shape  parameter  for  Gamma 
(   /^;  c.  j_,  /)  j_ )  prior  distribution 


A 

A- 


empirical  Bayes  estimate  of  scale  parameter  for  Gamma 
(   /.j_;  y.  ^,    f;^)    prior  distribution. 


It  is  of  interest  to  note  that  for  situations  in  which 
the  airlift  model  is  applicable,  problem  P2  is  just  a 
special  case  of  the  generalized  assignment  problem.  In  the 
formulation  of  the  generalized  assignment  problem 
[Ref.  24 :pp. 345-346] ,  each  agent  (or  man)  can  perform  or  be 
assigned  more  than  one  task  (or  job)  provided  that  the 
resource  available  to  the  agent  is  not  exceeded;  whereas,  in 
the  classical  assignment  problem  [Ref.  23:pp.  61-62],  each 
agent  can  be  assigned  to  only  one  task.  In  both 
formulations,  each  task  must  be  assigned  to  exactly  one  of 
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the  agents.  The  airlift  model  can  be  interpreted  as  a 
special  case  of  this  generalized  assignment  problem  in  which 
one  agent  must  perform  all  N  tasks,  however,  he  has  ( u^  +  1) 
alternative  ways  to  accomplish  each  task. 

Problem  P2  is  equivalent  to  PI  and  even  though  the 
number  of  variables  and  constraints  have  increased,  it  is 
the  preferred  formulation  because  it  eliminates  the  complex 
nonlinear  objective  function  of  PI.  The  coefficients  fj^^ 
can  be  directly  computed  using  recursive  formulas  for  the 
Poisson  or  negative  binomial  cumulative  distribution 
functions;  whereas,  Pj^(K^)  could  not  since  it  was  itself  a 
function  of  the  decision  variable.  Additionally,  because  of 
the  assumption  of  the  independence  of  failures,  other 
failure  distributions  can  be  readily  introduced  into  this 
linear  formulation.  Furthermore,  the  formulation  is  general 
enough  so  that  additional  constraints  can  be  incorporated  as 
required  and  as  will  be  seen  in  the  following  chapter,  it 
lends  itself  easily  to  programming  in  an  algebraic  modeling 
language. 
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IV.   SOLUTION  METHODOLOGY  AND  IMPLEMENTATION 

The  purpose  of  this  study  is  to  determine  an  optimal 
stockage  level  of  spare  secondary  reparables  for  each  of  the 
three  different  scenarios  described  in  the  introduction. 
Of  course  "optimality"  is  with  reference  to  the  failure 
model  adopted  and  the  relevance  of  the  data  used  to  estimate 
its  parameters.  The  determination  of  the  stockage  levels 
is  the  subject  of  this  chapter.  The  binary  linear  nature  of 
the  formulation  presented  in  the  preceding  chapter  suggests 
that  several  different  methods  of  solution  are  available  to 
determine  the  optimal  stockage  policy  for  a  given  scenario. 
Each  alternative  method  will  be  discussed  and  justification 
for  the  method  chosen  is  provided  in  addition  to  its 
implementation. 

A.   SOLUTION  APPROACHES  AND  METHODOLOGY 

The  mode  of  deployment  in  a  particular  scenario 
determines  the  number  of  constraints.  The  initial  solution 
approach  to  the  scenario  in  which  the  airlift  model  was 
applicable  utilized  the  technique  of  Lagrangian  relaxation. 
Although  it  is  standard  procedure  to  establish  bounds  for 
integer  programs  (such  as  problem  P2  )  using  a  linear 
programming  (LP)  relaxation,  it  was  initially  desired  to 
avoid  dependency  on  commercial  LP  solvers  not  locally 
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available  to  the  combat  service  support  planners  of  the 

Second  Force  Service  Support  Group  at  Camp  Lejeune,  North 

Carolina.   Therefore,  Lagrangian  relaxation  was  adopted 

because  it  lends  itself  nicely  to  programming  in  readily 

available  FORTRAN  77.    It  must  be  emphasized  that  the 

structure  of  problem  P2  did  not  demand  that  this  specific 

solution  technique  be  applied.   In  fact,  the  LP  relaxation 

of  problem  P2  could  have  been  solved  using  an  LP  solver. 

Thus,   Lagrangian  relaxation  was  utilized  not  out  of 

necessity,  but  primarily  because  it  offered  a  realistic 

alternative  to  the  use  of  a  commercial  LP  optimizer. 

Briefly  stated,  the  technique  of  Lagrangian  relaxation, 

as  outlined  by  Fisher  [Ref.  25:pp.  10-21],  involves  moving 

complicating  constraints  into  the  objective  function  using 

the  product  of  the  Lagrangian  multiplier  and  the  constraint 

violation  as  a  penalty  term.   In  particular,  examination  of 

problem  P2  in  the  context  of  the  airlift  scenario  results  in 

the   following  binary   linear   formulation   (actually   a 

generalized  assignment  formulation): 

INDICES: 

j -,  J    ^Y   secondary  reparables  or  components 

j  =  0,  ...,U;   number  of  secondary  reparables  i 
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DATA: 

fij  ln(Pi(j)) 


M- 


total  shipping  weight  for  j  spares  of  type  i 
[lbs} 


ACL  allowable  cargo  load  of  airplane  [lbs] 

W/        maximum   allowable   number   of   spares   for 
component  i 

DECISION  VARIABLE: 

_Y..  _  J^    if  j  spares  of  type  i  are  selected 
'■'        \0         otherwise 


(PB) 


FORMULATION: 

max 

s.t. 

/=ly=0 

H 

7=0 

,r^e{0,l} 
The  only  complicating  constraint  is  the  allowable  cargo  load 

(ACL)  constraint.   The  relaxed  formulation  with  sca.la.r 

Lagrangian  multiplier,  ).     ,    consists  of  N  separable,  multiple 

choice   problems  which  are  easy  to  solve.   On  the  other  hand, 

for  those  situations  in  which  both  the  weight  and.  volume 

constraints  of  problem  P2  must  be  considered   (i.e., 

Mediterranean  MAU  or  MPF  MAB  scenarios),  the  Lagrangian 

multiplier,  A  ,  is  now  vector-valued  (i.e.,  A  =(  X-^,    ).  2)    since 

two  constraints  must  be  relaxed,  and  this  is  no  longer  a 
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trivial  problem  to  solve.  The  primary  reason  for  this 
complication  is  that  the  solution  approach  now  involves  a 
two-dimensional  search  which  is  more  difficult  to  program 
than  the  simple  bisection  search  associated  with  the  single, 
scalar  Lagrangian  multiplier.  The  problem  is  not  impossible 
to  solve  and  methods  exist,  such  as  a  coordinate  search 
method  suggested  by  DeWolfe,  Stevens,  and  Wood  [Ref.26: 
pp. 8-10],  which  appears  to  offer  good  results  for  this  type 
of  problem. 

As  additional  constraints  (such  as  budget,  holding 
(storage)  cost,  or  shortage  ( stockout )  cost  constraints)  are 
added,  the  programming  difficulty  associated  with  the  search 
procedure  of  the  Lagrangian  function  becomes  prohibitive. 
Therefore,  to  accommodate  the  possible  inclusion  of 
additional  constraints  in  the  formulation  of  problem  P2,  it 
was  decided  to  utilize  the  Generalized  Algebraic  Modeling 
System  (GAMS)  in  conjunction  with  the  Zero/One  Optimization 
Methods  (ZOOM)  mixed  integer  program  solver.  The  use  of 
this  optimization  package  not  only  avoids  the  restrictive 
programming  considerations  discussed  above,  but  it  can  also 
incorporate  the  PEB  as  an  initial  step,  prior  to  solution  of 
P2  or  P3. 

The  approach  to  the  PEB  is  to  maximize  the  marginal 
likelihood  function  associated  with  each  secondary  reparable 
i  with  respect  to   a^^  and   /J^.   This  marginal  likelihood 
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function  is  a  nonlinear  function  of  the  two  decision 

variables.     The  formulation  for  this  maximization  problem 

is: 

INDICES: 
/■=!,...  ,iV  secondary  reparable  or  component 

m=l,...  ,A/j  source  of  data  (e.g.,  tank  serial  number  508282) 


DATA: 


im 


number  of  observations  of  component  i  failures 
for  data  source  m 


^im  total  observed  exposure  time  for  component  i  of 

data  source  m. 


«/ 


upper  bound  on  the  shape  parameter  of  the  Gamma 
prior  distribution  for  component  i. 


^i  upper  bound  on  the  scale  parameter  of  the  Gamma 

prior  distribution  for  component  i. 


DECISION  VARIABLES: 


a 


shape  parameter  for  Gamma  prior  distribution 
associated  with  compenent  i. 


^i  scale  parameter  for  Gamma  prior  distribution 

associated  with  component  i. 


FORMULATION: 


Ml 


A 


max 


m- 


im 


..  1  "J-^i^-d       ^    [,^  +  (i,    J      ^    l,^  +  P 


{PA) 


s.i.        0  <  a,  <  a,- 

Standard  maximization  by  taking  partial  derivatives  with 
respect  to    a^  and   /]j_  does  not  yield  a  satisfactory  closed 
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form  solution  which  can  be  practically  implemented.  Thus, 
this  nonlinear  maximization  problem  is  also  solved  using 
GAMS  in  conjunction  with  the  nonlinear  program  solver, 
MINOS  [Ref.  27].  These  values  are  then  incorporated  into 
the  optimization  of  problem  P2.  In  summary,  under  the 
current  solution  approach,  the  entire  model  can  be  coded  for 
computer  processing  on  either  a  mainframe  system  or  personal 
computer  and  solved  using  a  single  modeling  system  (i.e., 
GAMS)  and  its  compatible  solvers.  This  approach  is 
sufficiently  flexible  to  incorporate  additional  constraints 
or  decision  variables. 

B.   IMPLEMENTATION 

Since  the  decision  was  made  to  determine  the  optimal 
stockage  policy  for  all  scenarios  using  GAMS,  no  specifics 
of  the  original  solution  approach  involving  Lagrangian 
relaxation  will  be  mentioned  other  than  the  fact  that  it  is 
implemented  by  means  of  a  FORTRAN  77  computer  program.  The 
FORTRAN  program  is  based  upon  a  similar  one  developed  by 
DeWolfe  [Ref.  17:pp.  23-24].  A  copy  of  the  source  code 
developed  for  this  study  is  provided  in  Appendix  A  and 
examples  of  its  two  required  input  files  and  of  the  single 
output  file  are  given  in  Appendices  B  and  C,  respectively. 

The  solution  approach  adopted  centers  around  the  coding 
of  problems  P2  and  P4  for  computer  processing  using  the  GAMS 
language.   The  reader  is  referred  to  the  excellent  tutorial 
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by  Rosenthal  [Ref.  28]  for  the  specifics  of  this  versatile 
mathematical  modeling  language.  Upon  coding,  the  problems 
are  solved  using  the  GAMS  compatible  solvers  of  MINOS  and 
ZOOM  for  nonlinear  programs  and  mixed  integer  programs 
respectively.  Specifically,  the  source  code  provided  in 
Appendix  D  incorporates  the  successive  solve  feature  of  GAMS 
to  first  solve  the  nonlinear  maximization  problem 
represented  by  problem  P4.  The  results  of  this  maximization 
provide  the  estimates  of  the  shape  (  a^^)  and  scale  (  p^) 
parameters  of  the  prior  distributions  which  are  required 
before  the  second  binary  linear  maximization  (P2)  can  be 
performed.  The  final  result  is  the  optimal  stockage  level 
of  spare  secondary  reparables  for  the  given  scenario.  Even 
though  the  source  code  in  Appendix  D  is  specific  to  the 
airlift  model  ( P3 ) ,  GAMS  makes  it  very  easy  to  incorporate 
additional  constraints.  For  example,  if  a  volume  constraint 
is  also  applicable,  a  single  GAMS  PARAMETER  statement  to 
accomodate  the  volume  data  and  EQUATION  statement  to  reflect 
the  volume  constraint  are  the  only  additions  necessary. 

The  source  code  provided  is  compatible  with  either 
mainframe  or  personal  computer  usage.  When  run  on  the  IBM 
3033AP  mainframe  computer,  GAMS  version  2.05  was  utilized  in 
conjunction  with  MINOS  Version  5.1  and  ZOOM  version  2.1  and 
on  the  Zenith  Z-248  personal  computer  (equipped  with  80287-8 
math  coprocessor  and  640  Kbyte  memory  expansion  board),  PC 
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GAMS  Version  2.04  was  utilized  in  conjunction  with  MINOS 
Version  5.0  and  ZOOM  Version  2.1. 
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V.   RESULTS 

This  chapter  summarizes  results  obtained  using  actual 
operational  maintenance  data  to  determine  an  optimal  level 
of  spare  secondary  reparables  for  a  Marine  Amphibious  Unit 
(MAU)  in  a  realistic  exercise  scenario.  Verification  of  the 
model  is  accomplished  using  simulated  failure  data. 
Estimation  of  the  parameters  of  the  prior  distribution  by 
an  analysis  of  the  data  is  discussed  and  demonstrated;  this 
realizes  computationally  the  PEB  procedure.  Finally, 
computational  experience  for  both  the  mainframe  and  Zenith- 
248  personal  computer  is  briefly  discussed. 
A.   MODEL  TESTING 

The  data  used  in  the  development  and  testing  of  all 
models  was  provided  by  the  Operations  Section  of  the 
Supported  Activities  Supply  System  (SASSY)  Management  Unit 
in  coordination  with  the  Maintenance  Information  Systems 
Management  Office  of  the  Second  Force  Service  Support  Group 
at  Camp  LeJeune,  North  Carolina.  This  data  was  extracted 
from  the  ERO  History  File  of  the  Marine  Corps  Integrated 
Maintenance  Management  System  and  analyzed  according  to  the 
procedure  presented  in  Chapter  II.  This  analysis  provided 
the  basic  parameters  for  equation  (3.5)  (or  (3.4))  from 
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which  the  objective  function  coefficients  were  generated 
using  GAMS . 

The  airlift  model  ( P3 )  was  used  to  determine  the  optimal 
stockage  policy  for  twenty-five  different  secondary 
reparables  associated  with  the  pacing  items  of  a  Marine 
Amphibious  Unit  participating  in  a  realistic  short-term 
exercise.  This  problem  generates  190  binary  variables  and 
26  constraints.  The  recommended  stockage  level  for  the 
secondary  reparables  associated  with  the  three  pacing  items 
(i.e.,  M60A1  Tank,  LTP7A1  Amtrack,  and  MlOlAl  105mm 
Howitzer)  in  this  exercise  is  presented  in  Table  2.  Table  2 
displays  the  levels  computed  using  both  the  GAMS/  MINOS/ZOOM 
package   and   the   Lagrangian   relaxation   technique. 

For  this  exercise  scenario  99.89%  of  the  allowable  cargo 
load  was  utilized  and  the  objective  function  value  was 
within  0.07%  of  the  upper  bound  established  during  the 
branch  and  bound  procedure  used  by  the  ZOOM  solver.  The 
probability  of  mission  completion  given  the  stockage  level 
was  0.4149  and  this  highlights  one  of  the  advantages 
associated  with  using  this  model.  This  low  probability 
implies  that  the  allowable  cargo  load  must  be  increased.  The 
only  way  to  accomplish  this  is  to  increase  the  number  of 
cargo  planes  assigned  to  this  exercise.  This  provides 
Marine  Corps  planners  with  a  tangible  basis  on  which  to 
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TABLE  2 


MAU  AIRLIFT  SCENARIO  SPARE  STOCKAGE  LEVELS 


SPARE  STOCKAGE  LEVEL 


GAMS/ 

LAGRANGIAN 

NSN 

END  ITEM 

ZOOM 

RELAXATION 

2815001245387 

M60A1 

2 

2 

2520002241867 

M60A1 

3 

3 

2815002395819 

M60A1 

2 

2 

2520001549632 

M60A1 

2 

2 

2530000886657 

M60A1 

2 

3 

2520001184942 

M60A1 

2 

2 

2920010135802 

M60A1 

1 

1 

1010010703803 

M60A1 

1 

1 

1025019388105 

M60A1 

4 

5 

1010010720397 

M60A1 

1 

1 

1025010708993 

M60A1 

1 

1 

2815004303480 

LVTP7A1 

3 

4 

2520003973384 

LVTP7A1 

5 

5 

2520001443385 

LVTP7A1 

4 

4 

5805014593410 

LVTP7A1 

4 

5 

2520008949532 

LVTP7A1 

2 

2 

2910011714636 

LVTP7A1 

2 

3 

2530011509757 

LVTP7A1 

2 

2 

5865012207848 

LVTP7A1 

3 

3 

4320012035660 

LVTP7A1 

2 

2 

2920002317276 

LVTP7A1 

1 

1 

2530000886650 

LVTP7A1 

2 

2 

1005011855059 

MIOIAI 

5 

4 

1005011086434 

MIOIAI 

5 

5 

1005011457709 

MlOlAl 

5 

5 
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justify  their  requests  for  additional  airlift  assets  from 
the  Military  Airlift  Command  of  the  U.S.  Air  Force.  Table  3 
summarizes  the  results  for  both  solution  approaches. 


TABLE  3 


SUMMARY  OF  RESULTS  FOR  MAU  AIRLIFT  SCENARIO 


Solution 
Technique 


Objective 
Function 
Value 


Probability 
of  Mission 
Completion 


Percentage 
of  ACL  Used 


Deviation 
from 

Optimality 


GAMS/MINOS/ 
ZOOM 


-0.8797 


Lagrangian 
Relaxation    -0.8820 
( w/Heuristic ) 


0.4149 


0.4139 


99.89% 


100% 


0.07% 


0.47% 


B.   MODEL  VERIFICATION 

Verification  of  the  modeling  approach  was  conducted  by 
using  the  model  to  verify  the  fact  that  the  negative 
binomial  predictive  function  reduces  to  the  Poisson-Gamma 
function  as  the  number  of  failures  observed  (and  accordingly 
the  total  exposure  time)  approaches  infinity.  Thus,  the 
Poisson  function  provides  the  limiting  case  or  best  possible 
stockage  level  attainable  assuming  that  the  failure  rates 
for  each  secondary  reparable  are  precisely  known. 

To  accomplish  this  verification,  the  MAU  airlift  model 
was  run  using  the  Poisson  objective  function  (3.4)  with 
fixed  failure  rates.  The  resulting  stockage  levels 
represent  the  best  stockage  policy  that  can  be  achieved. 
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Then,  using  those  failure  rates,  twenty-five  different  sets 
of  simulated  failure  data  were  generated  for  each  of  10,  25, 
50,  100,  200  and  500  failures.  These  simulated  failure 
times  were  randomly  generated  using  the  probability  integral 
transform  to  obtain  exponential  failure  times  from  uniform 
(0,1)  deviates  provided  by  the  LLRANDOMII  Random  Number 
Generation  Package  [Ref.  29].  The  Bayesian  analysis  of  this 
data  was  performed  resulting  in  the  negative  binomial 
predictive  distribution  and  the  objective  function,  equation 
(3.5).  A  noninformative  Gamma  prior  distribution  with  a  =1 
and  /j  =1  was  employed  in  this  process.  Stockage  levels  for 
this  method  are  displayed  in  Table  4.  Note  that  about  100 
failure  observations  are  needed  before  the  stockage  levels 
approach  those  computed  using  the  actual  failure  rates 
(i.e.,  the  optimal  stockage  policy  given  perfect 
information).  Figure  1  is  a  plot  of  the  interquartile  range 
(IQR)  of  the  probabiltiy  of  mission  completion  associated 
with  each  failure  sample  size.  This  graph  demonstrates  that 
as  the  number  of  observed  failures  increases,  the 
variability  (represented  by  the  IQR)  in  the  probability  of 
mission  completion  decreases. 

In  addition  to  the  above  verification,  the  fact  that  the 
model  and  solution  approach  responded  predictably  to  a  wide 
range  of  different  failure  observations  and  exposure  times 
indicates  that  the  model  and  approach  tend  to  be  robust. 
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TABLE    4 


MAU   AIRLIFT    SCENARIO    SPARE    STOCKAGE    LEVELS 


SPARE    STOCKAGE    LEVEL 


NSN 


10 


NEGATIVE    BINOMIAL 
25      50         100         200      500 


POISSON 


Ui 

> 
0 
O 


2815001245387 
2520002241867 
2815002395819 
2520001549632 
2530000886657 
2520001184942 
2920010135802 
1010010703803 
1025019388105 
1010010720397 
1025010708993 
2815004303480 
2520003973384 
2520001443385 
5805014593410 
2520008949532 
2910011714636 
2530011509757 
5865012207848 
4320012035660 
2920002317276 
2530000886650 
1005011855059 
1005011086434 
1005011457709 


3  3        3  2  2  2 

4  4  4  3  3  3 
2           2        2          2  2  2 

2  2  2  3  3  3 
1112  2  2 
1112  2  1 
1111  11 
0          0        11  11 

4  4  4  4  4  4 
1111  11 
10        0           1  11 

5  4  4  4  4  4 
5           5        5           5  5  5 

3  4  4  4  4  4 
5  5  5  5  5  5 
2  2  2  2  2  2 
12        2           2  2           2 

2  2        2  2  2  2 

3  2        2  3  3  3 

3  2  2  2  2  2 
1111  11 
1111  12 

4  4        3  4  4  4 

5  5  5  5  5  5 
5           5        5           5  5  5 


2 
3 
2 
3 
2 
1 
1 
1 
4 
1 
1 
4 
5 
4 
5 
2 
2 
2 
3 
2 
1 
2 
4 
5 
5 
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Figure  1 . 


Variability  in  P( Mission  Completion)  with 
Number  of  Observed  Failures 
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C.   BAYESIAN  AND  EMPIRICAL  BAYES  ANALYSIS 

Bayesian  analysis  and  particularly  PEB  suggest 
themselves  when  the  failure  rates  of  different  units  (e.g., 
tank  engines  on  different  tanks  or  amtrack  transmissions  on 
different  amtracks )  are  sufficiently  alike,  but  individual 
experience  is  small,  so  that  "borrowing  strength"  by  a 
suitable  pooling  process  can  improve  individual  estimates. 
In  this  study,  the  failure  rate  associated  with  each 
secondary  reparable  (  /  j_ )  was  assumed  to  follow  a  Gamma 
prior  distribution.  This  prior  distribution  represents  the 
probability  that  the  unknown  value  of  the  failure  rate  lies 
in  various  regions  of  the  parameter  space  [Ref.  20:p.  260]. 
In  other  words,  incorporating  the  prior  distribution  enables 
a  tighter  bound  on  the  unknown  value  of  the  failure  rate  to 
be  obtained  than  would  be  possible  by  using  the  data 
pertaining  to  that  rate  alone.  The  primary  disadvantage  of 
the  ordinary  or  classical  Bayesian  approach  is  that  the 
prior  distribution's  parameters  (i.e.,  a^  and  /5^)  are 
assumed  to  be  known.  This  deficiency  is  overcome  in  this 
study  by  employing  a  PEB  approach,  in  which  the  entire  data 
set  (  (/?,>„,  tf^).  m=l,...,Mi  )  associated  with  each  secondary 
reparable  i  is  used  to  estimate  the  parameters  of  the  prior 
distribution.  These  estimates  are  then  used  in  the  standard 
Bayesian  analysis  to  compute  the  posterior  and  predictive 
distributions.     The  results  are  unavoidably  influenced  by 
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the  particular  choice  of  prior  used  (here  a  Gamma). 
Sentitivity  tests  have  not  been  conducted,  but  indication  is 
that  the  procedure  is  rather  insensitive. 

A  specific  indication  of  this  insensitivity  is  seen  by 
comparing  the  airlift  model  stockage  levels  obtained  when 
minimal  prior  information  is  included  to  the  levels  obtained 
when  the  empirical  Bayes  estimates  of  the  prior  distribution 
are  included.  The  results  of  this  comparison  are  summarized 
in  Table  5. 

The  levels  in  the  column  labeled  "NONINFORMATIVE  PRIOR" 
were  obtained  by  running  the  airlift  model  with  the  randomly 
generated  failure  times  (described  in  the  previous  section) 
associated  with  the  number  of  failures  sample  size  of  ten. 
A  "diffuse"  or  noninf ormative  prior  distribution 
(i.e.,  a  ^=  |]j_  =  l)  was  used.  A  diffuse  prior  distribution  is 
"informationless  in  a  relative  sense  and  means  only  that  it 
is  diffuse  relative  to  the  sample  information" . [Ref.  30:p. 
198]  The  same  airlift  model  was  then  run  using  empirical 
Bayes  estimates  for  the  prior  parameters  and  these  levels 
appear  in  the  column  labeled  "INFORMATIVE  PRIOR".  These 
levels   are  only  slightly  closer  to  the  Poisson  levels   than 
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TABLE  5 


EMPIRICAL  BAYES  ESTIMATION 


STOCKAGE  LEVELS 


u 
> 

0 

o 
< 


NEGATIVE 

BINOMIAL 

POISSON 

NONINFORMATIVE 

INFORMATIVE 

NSN 

PRIOR 

PRIOR 

2815001245387 

3 

2 

2 

2520002241857 

4 

4 

3 

2815002395819 

2 

2 

2 

2520001549632 

2 

2 

3 

2530000886657 

1 

1 

2 

2520001184942 

1 

1 

1 

2920010135802 

1 

1 

1 

1010010703803 

0 

0 

1 

1025019388105 

4 

4 

4 

1010010720397 

1 

1 

1 

1025010708993 

1 

1 

1 

2815004303480 

5 

5 

4 

2520003973384 

5 

5 

5 

2520001443385 

3 

4 

4 

5805014593410 

5 

5 

5 

2520008949532 

2 

2 

2 

2910011714636 

1 

1 

2 

2530011509757 

2 

2 

2 

5865012207848 

3 

3 

3 

4320012035660 

3 

3 

2 

2920002317276 

1 

1 

1 

2530000886650 

1 

1 

2 

1005011855059 

4 

4 

4 

1005011086434 

5 

5 

5 

1005011457709 

5 

5 

5 
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are  those  computed  when  minimal  prior  information  is 
considered. 

Even  though  these  stockage  levels  are  only  slightly 
closer  to  the  Poisson  levels  for  this  specific  instance,  by 
simply  using  the  available  data  for  10  failures  to  compute 
empirical  Bayes  estimates  of  the  prior  distribution 
parameters,  spare  stockage  levels  comparable  to  those 
obtained  in  the  previous  section  using  failure  sample  sizes 
of  between  50  to  100  can  be  achieved.  This  is  particularly 
significant  since  we  have  no  control  over  the  number  of 
failures  observed  for  each  secondary  reparable  or  their 
exposure  times  and  would  not  usually  observe  this  number  of 
failures  (i.e.,  50-100)  during  the  reporting  period  covered 
by  the  ERG  History  File.  Also,  since  operational 
reliability  data  may  not  be  available  for  the  exact  location 
where  the  MAGTF  is  deploying,  pooling  the  data  from  sources 
which  have  operated  in  similar  environments  has  the 
advantage  of  providing  more  realistic  failures  rates  than  if 
this  prior  information  were  not  considered. 

D.   COMPUTATIONAL  EXPERIENCE 

So  as  to  keep  the  implementation  and  solution  of  the 
models  developed  in  this  study  as  general  and  flexible  as 
possible  existing  software  packages  (i.e.,  STATGRAPHICS, 
GAMS,  etc. )  were  utilized.  The  goal  was  to  obtain  solutions 
in  a  relatively  user-friendly  personal  computer  environment 
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requiring  a  minimum  of  operator  interaction.  Therefore, 
computational  efficiency  in  terms  of  computer  run  time  (or 
algorithmic  complexity)  was  not  of  primary  concern  and  is 
included  here  solely  for  completeness.  No  attempt  was  made 
to  employ  special  algorithms  that  exploit  problem  structure. 
This  is  not  to  downplay  the  importance  of  this  factor,  but 
ir  rather  to  emphasize  that  during  this  study  the  primary 

^,  concern  was  to  obtain  solutions  to  real-world  problems  as 

u. 

easily  and  efficiently  as  possible  with  a  minimum  of 
u. 

2  programming  effort. 

c: 

hi  Running  time  for  the  airlift  model  on  the  IBM  3033AP 

o 

o  using  the  GAMS/MINOS/ZOOM  package  was  3.9  seconds  in 

< 

ui 
u 

D 

a 

0 

a. 

hi 


contrast  to  the  176   seconds  obtained  when  run  on  the  Zenith 
Z-248  PC.   It  is  possible  on  the  IBM  3033AP  to  combine  the 
PEB  calculations  with  the  maximization  of  problem  P3  by 
c  using  the  successive  solve  feature  of  GAMS;  however,  this 

cannot  be  done  on  the  Zenith  Z-248  PC  since  internal 
parameters  in  the  PC  version  of  GAMS  are  exceeded. 
Therefore,  each  of  the  empirical  Bayes  maximizations  must  be 
performed  separately  when  using  this  PC.  The  resulting 
estimates  have  to  be  manually  input  into  the  appropriate 
GAMS  PARAMETER  statements  before  maximization  of  problem  P3 
can   be   accomplished. 

Even  though  the  Lagrangian  relaxation  procedure  with 
heuristic  was  not  adopted  as  the  solution  technique   for  the 
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reasons  discussed  in  Chapter  IV,  it  provides  much  faster 
solution  times  to  the  airlift  model  ( P3 )  on  both  the  IBM 
3033AP  and  Zenith  Z-248  PC.  The  Lagrangian  procedure,  coded 
in  FORTRAN  77,  provides  a  running  time  of  0.32  seconds  on 
the  IBM  3033AP  and  about  30  seconds  on  the  Zenith  Z-248  PC 
where  it  was  coded  using  Microsoft  FORTRAN  (Version  3.2). 

Finally,  comparison  of  the  stockage  levels  obtained 
using  GAMS/MINOS/ZOOM  with  those  obtained  using  LP 
relaxation  for  each  of  the  150  simulated  failure  time  cases 
revealed  that  in  every  case  the  LP  relaxation  provided 
identical  stockage  levels  with  the  exception  that  it 
fractionated  for  exactly  one  of  the  twenty- five  secondary 
reparables.  So,  if  a  general  policy  was  adopted  to  round 
down  to  the  lower  of  the  fractionated  stockage  levels  (so  as 
not  to  violate  the  ACL  contrast),  the  LP  relaxation  would 
provide  comparable  results  to  those  obtained  using  the  mixed 
integer  program  solver,  ZOOM.  This  is  an  important  result 
because  the  run  times  associated  with  the  LP  relaxation  were 
significantly  less  than  those  using  ZOOM.  In  particular, 
if  the  rounding  down  policy  is  applied  to  the  LP  relaxation 
of  the  airlift  model  in  this  study,  the  resulting  integer 
stockage  level  provides  an  objective  function  value  which  is 
within  1.87%  of  the  LP  upper  bound  or  1.73%  of  the  upper 
bound  established  by  the  branch  and  bound  procedure  utilized 
by  ZOOM.   This  is  not  quite  as  good  as  the  result  obtained 
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directly  from  ZOOM  (within  0.07%  of  optimality ) ,  however, 
for  those  situations  in  which  computer  run  time  may  be  a 
consideration  (e.g.,  as  problem  size  increases)  or  in  which 
very  tight  bounds  may  not  be  required,  a  very  good  solution 
can  be  obtained  from  the  much  simpler  and  less  costly  LP 
relaxation. 
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VI.   CONCLUSIONS  AND  RECOMMENDATIONS 

In  this  thesis,  the  problem  of  determining  spare 
stockage  levels  for  deployable  MAGTFs  was 
solved  by 

(a)  stochastically  modeling  demand  utilizing  operational 
data  from  the  Marine  Corps  Integrated  Maintenance 
Management  System,  and  then 

(b)  using  mathematical  programming  to  optimize  a 
suitable,  operationally  relevant  objective  function. 

To  accomplish  this,  a  new  measure  of  effectiveness  was 

defined  which  provides  the  operational  commander  on  the 

ground  with  an  objective  indicator  of  the  overall  impact  of 

the   stockage   policy   on   the   probability   of   mission 

completion.   This  measure  of  effectiveness  was  evaluated  as 

the   end  result  of   a  probabilistic  modeling   step,   a 

statistical  data  analysis,  and  a  subsequent  optimization. 

The  generality  of  the  model  formulation   makes  it 

adaptable  to  any  of  the  three  scenarios  of  present  interest 

to  the  Marine  Corps.   The  use  of  a  PC-compatible  algebraic 

modeling  language  (GAMS)  makes  it  particularly  easy  to  add 

or  delete  constraints  as  required.   Specifically,  the  model 

and  solution  approach  were  applied  to  a  realistic  exercise 

scenario  and  an  essentially  optimal  stockage  policy  was 

determined  which  maximized  the  probability  of  mission 
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completion  within  the  constraints  imposed  by  the  deployment 
mode. 

During  the  analysis  of  the  MIMMS  data  base  (i.e.,ERO 
History  File)  numerous  problems  and  inconsistencies  were 
encountered.  As  a  result,  the  estimates  of  the  failure 
rates  for  some  of  the  secondary  reparables  are  higher  than 
would  be  the  case  if  all  the  potential  times  to  failure 
could  be  extracted  from  the  data  base.  To  overcome  this 
discrepancy,  it  is  highly  recommended  that  the  quarterly 
update  procedure  described  in  Chapter  II  be  conducted  at  the 
Force  Service  Support  Group  level.  If  this  cannot  be  done, 
then  at  the  very  least  the  meter  readings  associated  with 
all  pacing  items  should  be  recorded  and  input  into  MIMMS 
prior  to  each  quarterly  update.  This  would  insure  that 
beginning  and  ending  meter  readings  are  available  for  each 
serialized  pacing  item  during  the  reporting  period  covered 
by  the  ERO  History  File  and  this  would  prevent  the  loss  of 
valuable  exposure  time  data. 

The  limited  strategic  lift  assets  available  to  deploying 
Marine  Corps  units  demands  that  effective  utilization  be 
made  of  these  resources.  Since  only  a  small  percentage  of 
these  assets  can  be  allocated  for  spare  parts  inventories, 
emphasis  must  be  placed  on  stocking  those  parts  which  would 
have  the  greatest  impact  on  the  completion  of  the  assigned 
mission.   The  reliability-based  model  developed  in  this 
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study  accomplishes  this  objective  by  maximizing  the 
probability  of  mission  completion  within  the  constraints  of 
the   mode   of   deployment. 
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APPENDIX    A 
LISTING    OF    SOURCE    CODE    FOR    LANGRANGIAN    RELAXATION 

PROGRAM    RELAX2 

CHARACTER      N5NX13, N0MENX7 

INTEGER      NC0MP,NEQUIP,NSPARE,MAXSPR,IDNUM,BSPARE,SPBND,NDEN5, 
XNaPER,NXPOSE 

REAL5(8         P,N,MAXWT,DINF,OPTDEV,ZLOW,ZUP,PSUC,HILAMD,WT,MTBF, 
'"^  XDURMIS,ZHLON,XPOSE, ALPHA, BETA 

t  PARAMETER    (  NC0MP  =  25,  NEQUIP  =  5,  MAXSPR  =  1  0  ,  MAXWT  =  ^'+^00  .  ODO, 

-'  XDINF  =  1.0D20) 

DIMEfJSION      NSN(NC0MP),N0MEN(NCOMP),NSPARE(NC0MP),IDNUM(NCOMP), 
-  XP(fJCOMP,0:  MAXSPR),K'CfJCCMP,0:MAXSPR),  B5PARE(  f.'COMP  )  ,  SPBNDC  NCOMP) , 

ir  5(:JDENS(NC0MP),N0PERCfJC0MP),XP0SE(NC0MP),NXP0SE(NC0MP), 

^  5tDURMIS(NC0MP),ALPHACNC0MP),BETA(NC0MP) 

S  CALL  INIT2(DINF,MAXSPR,NEQUIP,MAXWT,P,W,HSN,N0MEN,HILAMD.NSPARE, 

>  XSPBND, NCOMP, I DNUM,NDENS, NOP ER,XP0SE,NXP05E,DURMIS, ALPHA, BETA) 

0 

^  CALL  BGUND(DINF,P,W,NCOMP,MAXSPR, SPEND, MAXWT, HI L AMD, NSPARE,ZLOW, 

I-  XZUP,WT,B5PARE) 

< 

Q  CALL  HERIST(DINF,P,W,NC0MP,MAXSPR,ZHL0W,WT,MAXWT,SPBND,3SPARE, 

u  J^LARRAY) 

u 

Q 
O 

c: 

^  OPTDEV  =  -l.ODO  X  100. ODO  x  (l.ODO  -  (ZHLOW/ZUP)) 

c: 


WRITE  (10,50) 
50    FORMAT  (  'I'/'O',  •  — 


WRITE  (10,100)  ZHLOH 
100   FORMAT  (IX,  'OBJECTIVE  FUNCTION  VALUE  IS  :',D17.10) 

WRITE  (10,150)  OPTDEV 
150   FORMAT  CO'. 'THIS  SOLUTION  IS  WITHIN  ',D17.10,'  '/.    OF  THE  OPTIMAL') 

PRINT  X, 'ZUP=',ZUP, 'ZLOW=',ZLON, ' ZHLOW= ' , ZHLOW 

P5UC  =  EXP(ZHLOW) 

WRITE  (10,175)  PSUC 
175   FORMAT  ( ' 0 ' , • P(MISSION  SUCCESS)  =  •,F7.5) 

WRITE  (10,180) 

130   FORMAT  (  '0',  ' S 

X' '  ) 

WRITE  (10.200) 
200   FORMAT  ( ' 0 ' , 5X, 'NSN' , lOX, • END  ITEM' , 6X, ' SPARE  LEVEL') 

DO    250    I  =  l,riCOMP 

WRITE  (10,300)  NSN(I),NOMEN(IDNUM(I)),BSPARE(I) 
250   CONTINUE 
300   FORMAT  ( ' 0 ' , A13, 7X, A7 , 7X, 14) 

PRINT  X, 'WEIGHT  USED= ' , WT, 'MAXWT= ' , MAXWT 

STOP 
END 
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SUBROUTINE  INITZC DINF, MAXSPR, NEQUIP, MAXWT, P . W, NSN, NOMEN , HILAMD, 
5(N3PARE,  SPEND,  NCOMP,  I DNUM,NDENS,N0PER,XP0SE,NXPC5E,DURMIS,  ALPHA, 
XBETA) 

INTEGER   NCOMP 
REALMS    MAXWT 

INTEGER   I DNUM( NCOMP ),NDENS( NCOMP ),NOPER( NCOMP ),LMT,MAXSPR, 
^NSPARE(NC0MP;,SPBNDCNC0MPJ,NEQUIP,NXPQ5ECNC0MP) 

REAL?€S    W(  NCOMP,  0:MAXSPR),DURMIS(  NCOMP)  ,  DINF,  PI ,  Q,  BETA  (  NCOMP  ) , 
*PT  EMP1,P(  NCOMP,  O-.MAXSPR),  HI  L  AMD,  HLTEMP,CUMSUM,XPOSE(  NCOMP  ),PP, 
XALPHA(NCOMP) 

CHARACTER   NSNX13( NCOMP) , N0MENX7 ( NCOMP) 

READ  (1,10)  (NSN(I),IDNUM(I),XP0SE(I),NXP0SE(I),WCI,1),SPBND(I), 
XALPHACI),BETA(I),  1=1, NCOMP) 
10    FORMAT  (A13,2X,I2,2X,F9.1,2X,I4,2X,F6.1,2X,I3,2X,F7.4,2X,F8.2) 

READ  (2,20)  (NOMEN(I),NDENS(I),NOPER(I),DURMIS(I),  I=1,NEQUIP) 

20  FORMAT  ( A7 , ^X, 13, ^X, 13, ^X, F6 . 1 ) 

DO  21  1=1, NCOMP 
PRINT  )f,NSN(I),IDNUM(I),XPOSE(I),NXPOSE(I),W(I,l),SPBr<D(I), 
?eALPHACI),BETA(I) 

21  CONTINUE 

DO  22  I=1,NEQUIP 
PRINT  ?(,NOMEN(I),NDENS(I),NOPER(I),DURMIS(I) 

22  CONTINUE 

DO  25  1=1, NCOMP 
NSPARE(I)  =  0 
DO  50  J=0,MAX5PR 
P(I,J)  =  DINF 

IF  (J  .NE.  1)  W(I,J)  =  DINF 
50       CONTINUE 
25    CONTINUE 

DO  75  1=1, NCOMP 

PI  =  (XP05ECI)+BETA(I))/(NDENS(IDNUM(I))X 
XDURMISdDNUMd))  +  XPOSE(I)  +  BETACD) 
Q  =  1 .ODO  -  PI 

PP  =  PI  J^x  (NXPOSE(I)  +  ALPHA(I)) 
LMT  =  NDENSdDNUMd))  -  NOPER(IDNUM(  I ) ) 
DC  100  J  =  0,SPBNDd) 

IF  (J  .NE.  1)  Wd,J)  =  JXW(I,1) 
CUM3UM  =  O.ODO 
DO  125  K=0,(LMT  +  J) 
IF  (K  .EQ.  0)  THEN 

PTEMPl  =  l.ODO 
ELSE 
PTEMPl  =  PTEMPl)(((((NXPOSEd)  +  ALPHA(I))  +  (K-l);xi.0D0)/K)xQ 
END  IF 

CUMSUM  =  CUMSUM  +  PTEMPl 
125         CONTINUE 

CUMSUM  =  CUMSUM  x  PR 
Pd,J)  =  LOG(CUMSUM) 

PRINT  X,'I=',I,'J=',J,'LN=',P(I,J), •CUMSUM=', CUMSUM 
100      CONTINUE 
75    CONTINUE 
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HILAMD  =  -DINF 
DO  200  I=1,NC0MP 

DO  225  J=1,SPBND(I) 

HLTEMP  =  (  P(I,J)  -  P(I,0)  )  /  W(I,J) 
IF  (HLTEMP  .GT.  HILAMD)  HILAMD  =  HLTEMP 
225      CONTINUE 
200   CONTINUE 

PRINT  X, •HILAMBDA=', HILAMD 

RETURN 

END 

SUBROUTINE  BOUND( DINF, P, W, NCOMP, MAXSPR, SPEND, MAXWT, HILAMD, NSPARE, 
XZL01s,ZUP,NT,3SPARE) 

INTEGER   NCOMP,BSPARE(NCOMP),MAXSPR, SPEND (NCOMP),NSPARE(NCOMP) 

REAL3t8         BSTLAM,ZLOW.WT,P(NCOMP,0:MAXSPR)  .  W(  NCOMP,  0  :  MAXSPR)  ,  DINF, 
5(MAXWT,HILAMD,EPS,LFTEND,RTEND,2LBEST,LAMBDA,ZUP 

EPS  =  HILAMD/1.0D04 
LFTEND  =  O.ODO 
"  RTEND   =  1  .IDO^HILAMD 
ZLBEST  =  -DINF 
50    LAMBDA  =  (LFTEND  +  RTEND)/2.0D0 

CALL  MAXFCTC DINF, P,W,NCOMP,MAXSPR, SPEND, MAXWT, LAMBDA,  NSPARE,ZLOW, 
XZUP,WT) 

IF  (WT  .LE.  MAXWT)  THEN 
RTEND  =  LAMBDA 
IF  CZLOW  .GE.  ZLBEST)  THEN 
ZLBEST  =  ZLOW 
BSTLAM  =  LAMBDA 
END  IF 
ELSE 

LFTEND  =  LAMBDA 
END  IF 

IF  (  (RTEND  -  LFTEND)  . GT .  EPS  )  GOTO  50 

CALL  MAXFEA( DINF, P,W,NCOMP,MAXSPR, SPEND, BSTLAM, BSPARE,WT, ZLOW) 

RETURN 
END 
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SUBROUTINE  MAXFCTC DINF, P, W, NCOMP,MAXSPR, SPBND, MAXWT, LAMBDA, NSPARE, 
XZL0W,ZUP,WT) 

INTEGER   NCOMP,NSPARE(NCOMP),MAXSPR,SPBND(NCOMP), INDEX 

REALMS    LAMBDA, ZLOW,WT,P(NCOMP,0:MAXSPR),W(NCOMP,0:MAXSPR),DINF, 
XCT0TAL,CMPMAX,ZVALJ,0BJFCT,0BJVAL,ZUP,MAXWT 

WT  =  O.ODO 
CTOTAL  =  O.ODO 
ZLOW  =  O.ODO 

DO  25  I=1,NC0MP 
CMPMAX  =  -DINF 
INDEX  =  0 

DO  50  J=0,SPBNDCI) 

ZVALJ   =  P(I,J)  -  LAMBDA3(W(I,  J) 
OBJVAL  =  PCI, J) 
IF  (ZVALJ  .GE.  CMPMAX)  THEN 
NSPARE(I)  =  J 
CMPMAX  =  ZVALJ 
OBJFCT  =  OBJVAL 
INDEX   =  J 
END  IF 
50       CONTINUE 

CTOTAL  =  CTOTAL  +  CMPMAX 
ZLQH  =  ZLOW  +  OBJFCT 
WT  =  WT  +  W( I, INDEX) 
25    CONTINUE 

ZUP  =  CTOTAL  +  LAMBDA^MAXWT 

RETURN 
END 

SUBROUTINE  MAXFEAC DINF, P, W, NCOMP, MAXSPR, SPEND, B3TL AM, BSPARE, WT, 
XZLOW) 

INTEGER   NCOMP,BSPARE(NCOMP),MAXSPR,SPBND(NCOMP), INDEX 

REALX8    BSTLAM,ZLOW,WT,P(NCOMP,0:MAXSPR),W(NCOMP,0:MAXSPR),DINF, 
XCMPMAX,ZVALJ, OBJFCT, OBJVAL 

WT  =  O.ODO 
ZLOW  =  O.ODO 

DO  25  I=1,NC0MP 
CMPMAX  =  -DINF 
INDEX  =  0 

DO  50  J=0,SPBND(I) 

ZVALJ   =  P(I,J)  -  BSTLAM)(W(I,  J) 
OBJVAL  =  P(I,J) 
IF  (ZVALJ  .GE.  CMPMAX)  THEN 
BSPARE(I)  =  J 
CMPMAX  =  ZVALJ 
OBJFCT  =  OBJVAL 
INDEX   =  J 
END  IF 
50       CONTINUE 

ZLOW  =  ZLOW  +  OBJFCT 
WT  =  WT  +  W( I, INDEX) 
25    CONTINUE 

RETURN 
END 
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SUBROUTINE  HERISTC DINF, P, W, NCOMP, MAXSPR, ZHLOW, WT, MAXWT, SPBND, 
XBSPARE.LARRAY) 

INTEGER   NCOMP,BSPARE(NCOMP),MAXSPR,SPBND(NCOMP), INDEX 

REALX8    DINF,P(NCOMP,0:MAXSPR),W(NCOMP,0:MAXSPR),WT1,WT,MAXWT, 
XXSSWT,CHKVAL,L ARRAY (NCOMP),NUMER,DENOM,L BEST, ZHLOW 

CHKVAL  =  DINF/l.lDO 
HTl  =  WT 
25    XSSWT  =  MAXWT  -  WTl 

DO  50  I=1,NC0MP 

IF  ((BSPARE(I)  .GE.  SPBND(I))  .OR.  (P( I , ( BSPAREC I)+l) )  . GT . 
XCHKVAD)  THEN 

LARRAY(I)  =  -DINF 
GOTO  50 
END  IF 

NUMER  =  P(I,(BSPARE(I)+1))  -  P( I , BSPAREC I) ) 
DENOM  =  WCI, (BSPARE(I)+1))  -  W( I , BSPAREC I) ) 
IF  ((DENOM  .GT.  XSSWT)  .OR.  (NUMER  .LE.  l.OD-10))  THEN 
LARRAY(I)  =  -DINF 
GOTO  50 
ELSE 

LARRAY(I)  =  NUMER/DENOM 
END  IF 

IF  (LARRAY(I)  .LT.  O.ODO)  LARRAY(I)  =  -LARRAY(I) 
50    CONTINUE 

LBEST  =  -DINF 
DO  75  I=1,NC0MP 

IF  (LARRAY(I)  .GT.  LBEST)  THEN 
LBEST  =  LARRAYCI) 
INDEX  =  I 
END  IF 
75    CONTINUE 

IF  (LBEST  .LT.  (-CHKVAL))  GOTO  100 

WTl  =  WTl  -  W(INDEX,BSPARE(INDEX))  +  W( INDEX, ( BSPAREC INDEX)+1 ) ) 
BSPARECINDEX)  =  BSPAREC  IfJDEX)  +  1 
IF  (WTl  .LE.  MAXWT)  GOTO  25 
100   ZHLOW  =  O.ODO 
WT  =  O.ODO 
DO  125  I=1,NC0MP 

ZHLOW  =  ZHLOW  +  P( I , BSPARE( I ) ) 
WT  =  WT  +   W(I,BSPARE(I)) 
125   CONTINUE 

RETURN 
END 
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APPENDIX  B 


SAMPLE  INPUT  FILES 


PRINCIPAL  END  ITEM  FILE 


M60A1 

5 

^ 

60.0 

LVTP7A1 

12 

11 

40.0 

MIOIAI 

5 

4 

2000.0 

TABLE  B-1 


PRINCIPAL  END  ITEM  FILE  SPECIFICATIONS 


LINE 


All 


COLUMN 

!-• 

7 

8-: 

LI 

12- 

-14 

15- 

-18 

19- 

-21 

22- 

-25 

26- 

-31 

DATA  DESCRIPTION 


DATA  TYPE 


Nomenclature  of  end  item    Character 


Blank 

Quantity  of  end  item 

Blank 

Number  of  end  items 
required  operational 

Blank 

Mission  duration 


Integer 


Integer 


Real*8 
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SECONDARY  REPARABLE  COMPONENT  FILE 


28150012<45387 

2091 

0 

7 

2000 

0 

5 

0 

34000 

58 

6903 

2520002241867 

2091 

0 

8 

1000 

0 

5 

1 

00000 

1 

0000 

2815002395819 

2091 

0 

3 

500 

0 

5 

1 

00000 

1 

0000 

2520001549652 

2091 

0 

5 

500 

0 

5 

0 

10110 

13 

1350 

2530000836657 

2091 

0 

1 

100 

0 

10 

1 

00000 

1 

0000 

2520001184942 

2091 

0 

1 

250 

0 

5 

1 

ocooo 

1 

0000 

2920010135802 

2091 

0 

1 

500 

0 

10 

1 

00000 

1 

0000 

1010010703305 

2091 

0 

1 

1000 

0 

10 

1 

00000 

1 

0000 

1025019388105 

2091 

.0 

5 

100 

0 

5 

1 

00000 

1 

0000 

1010010720397 

2091 

.0 

3 

2000 

0 

5 

1 

00000 

T^ 

0000 

102501070S993 

2091 

0 

1 

1000 

0 

5 

1 

i 

00000 

i 

0000 

2S15004303480 

4. 

3253 

0 

8 

500 

0 

5 

0 

44330 

142 

9404 

2520003973384 

2 

3258 

0 

23 

1000 

0 

5 

0 

69700 

69 

0593 

2520001443385 

2 

3253 

0 

14 

2000 

0 

5 

1 

00000 

1 

0000 

5805014593410 

2 

3253 

0 

10 

250 

0 

5 

0 

20510 

27 

7348 

2520008949532 

2 

3258 

.0 

3 

500 

0 

5 

1 

00000 

1 

0000 

2910011714636 

2 

3253 

.0 

1 

100 

0 

10 

1 

00000 

1 

0000 

2530011509757 

2 

3258 

.0 

4 

1000 

0 

10 

1 

.00000 

1 

0000 

5865012207848 

2 

3258 

0 

2 

100 

0 

10 

1 

ocooo 

1 

0000 

4320012035660 

2 

5253 

0 

3 

500 

0 

10 

1 

00000 

1 

.0000 

2920002317276 

2 

3258 

0 

1 

500 

0 

10 

1 

00000 

1 

.0000 

2530000886650 

2 

3258 

0 

1 

250 

0 

5 

1 

.00000 

1 

.0000 

1005011355059 

3 

20996 

0 

3 

1000 

0 

5 

1 

00000 

1 

0000 

1005011086434 

3 

20996 

0 

11 

250 

0 

5 

1 

00000 

1 

0000 

1005011457709 

3 

20996 

0 

7 

500 

0 

5 

1^ 

00000 

1 

0000 
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TABLE  B-2 

SECONDARY  REPARABLE/COMPONENT  FILE  SPECIFICATIONS 

LINE      COLUMN     DATA  DESCRIPTION  DATA  TYPE 

All       1-13      National  Stock  Number  (NSN)    Character 

14-15      Blank 

16-17      Identification  Number         Integer 
(1  for  M60A1,  2  for  LVTP7A1, 
3  for  MIOIAI) 

18-19      Blank 

20-28      Total  exposure  time  for       Real*8 
component 

29-30      Blank 

31-34     Number  of  observed  failures    Integer 
for  component 

35-36      Blank 

37-42      Component  shipping  weight      Real*8 

43-44      Blank 

45-47     Maximum  possible  stockage      Integer 
for  component 

48-49      Blank 

50-56      Shape  parameter  of  prior      Real*8 
distribution  for  component 
failure  rate 

57-58      Blank 

59-66      Scale  parameter  of  prior      Real*8 
distribution  for  component 
failure  rate 
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APPENDIX    C 

SAMPLE    OUTPUT    FILE 
OBJECTIVE    FUNCTION    VALUE    IS    •• -0  .  8S2026315SD+00 
THIS    SOLUTION    IS   WITHIN      0 .  4761377567D+00    '/.   OF   THE   OPTIMAL 
P(MISSION    SUCCESS)    =    0.^41394 


NSN 

r 

2S150012453S7 
~  2520002241S67 

2815002395819 

ir  25200015<+9632 

\  2530000S86657 

>  252000118^+942 

O 

O  2920010135802 

<  1010010703803 

2  1025019322105 

u 

1010010720397 

1025010708993 

2815004303480 

2520003973334 

2520001443385 

5805014593410 

2520008949532 

2910011714636 

2530011509757 

5365012207848 

4320012035660 

2920002317276 

2530000886650 

1005011855059 

1005011086434 

1005011457709 


END  ITEM 

SPARE  LEVEL 

M60A1 

2 

M60A1 

3 

M60A1 

2 

M60A1 

2 

M60A1 

3 

M60A1 

2 

M60A1 

1 

M60A1 

1 

M60A1 

5 

M60A1 

1 

M60A1 

1 

LVTP7A1 

4 

LVTP7A1 

5 

LVTP7A1 

4 

LVTP7A1 

5 

LVTP7A1 

2 

LVTP7A1 

3 

LVTP7A1 

2 

LVTP7A1 

3 

LVTP7A1 

2 

LVTP7A1 

1 

LVTP7A1 

2 

MIOIAI 

4 

MIOIAI 

5 

MlOlAl 

5 
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APPENDIX  D 
GAMS  SOURCE  CODE 


SETS 

I   end  item  /E1875, E03^6 , E06^0/ 

J   component  /15^25/ 
K   spare  stockage  level   /SO^^SIO/ 

L   dummy  index  /DO^Dll/; 

SCALARS 

WTCAP   total  weight  capacity  (lbs)  /'4<4^00/ 
PSUC    probability  of  mission  completion; 

PARAMETERS 

EQDENS(I)  density  of  end  item  i 

/E1875  5 

E0S';6  12 

E06^0  5/; 

PARAMETERS 

NOPER(I)  number  end  item  i  required  operational  during  mission 

/E1S75  ^ 

E08A6  11 

E06^0  4/; 

PARAMETERS 

DURMIS(I)  mission  duration  for  end  item  i  (hours  rounds  miles) 

/E1S75  60 

E0846  ^0 

E0640  2000/; 

PARAMETERS 

SHIPWT(J)  shipping  weight  of  component  j    (lbs) 

/I  2000 

2  1000 

3  500 

4  500 

5  100 

6  250 

7  500 

8  1000 

9  100 

10  2000 

11  1000 

12  500 

13  1000 
1<;  2000 

15  250 

16  500 

17  100 

18  1000 

19  100 

20  500 

21  500 

22  250 

23  1000 

24  250 

25  500/; 
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PARAMETERS 

SPRMAX(J) 

max  number  of  spares  allowed 

/I 

5 

2 

5 

3 

5 

4 

5 

5 

10 

6 

5 

7 

10 

8 

10 

9 

5 

10 

5 

11 

5 

12 

5 

13 

5 

1"^ 

5 

15 

5 

16 

5 

17 

10 

13 

10 

19 

10 

20 

10 

21 

10 

22 

5 

23 

5 

24 

5 

25 

5/; 

PARAMETERS 

NOBS(J) 

number  of  data  points  observed 

/I 

7 

2 

8 

3 

3 

4 

5 

5 

1 

6 

1 

7 

1 

8 

1 

9 

5 

10 

3 

11 

1 

12 

8 

13 

23 

14 

14 

15 

10 

16 

3 

17 

1 

18 

4 

19 

2 

20 

3 

21 

1 

22 

1 

23 

3 

24 

11 

25 

7/; 

for  component  j 


<:. 
Z 

u 

>  NOBS(J)   number  of  data  points  observed  for  component  j 

O 

O 

< 

Q 
UJ 

u 

D 
Q 
O 

ir 

a. 
lii 
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PARAMETERS 

ALPHA(J) 

/I 

2 

3 

^  0 
5 
6 
7 
8 
9 

10 
11 

12  0 

13  0 
14 

15     0 

16 

17 

15 

19 

20 

21 

22 

23 

2'4 

25 


shape 

Q.3<i 


paramtr  for  Gamma  prior  distribution 


1011 


4433 
6970 

2031 


/ 


PARAMETERS 
BETA(J) 
/I 
2 
3 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 


scale  paramtr  for  Gamma  prior  distribution 
58.6903 


.1850 


14 
6 


.9404 

.0598 

.7348 


/; 


81 


PARAMETERS 

OK(I,J)   equals  1  if  component  j  is  a  component  of  end  item  i 
/E1875. (15(11)      1 
E0S46  .(12)(22)     1 
E06'40.  (23^25)     1/; 

PARAMETERS 

TXPOSE(I)  total  observed  exposure  for  end  item  i  (hrs  rds  mi) 

/E1S75  2091 

E08^6  3258 

E0640  20996/; 

PARAMETER 

LMT(I,J)   max  number  of  failures  for  component  j; 
LMT(I,J)$(  OK(I,J)  EQ  1  )  =  EQDENS(I)  -  NOPER(I); 

SETS 

OKLVL 1 ( I , J , K)   feasible  spare  stockage  levels  for  component  j; 
OKLVLKI,  J,K)$((0RD(K)-1  LE  SPRMAXCJ))  AND  (OKCI,J)  EQ  D)  =  YES; 

X   The  following  set  is  necessary  to  obtain  the  cumulative  distribution 
^      function  of  the  negative  binomial  density. 

SETS 

0KLVL2(I, J,K,L); 

QKLVL2(I, J,K,L)$(  C0RD(K)-1  LE  SPRMAX(J))  AND  (OK(I,J)EQ  1)  AND 
(ORD(L)-l  LE  LMT(I,J)+0RD(K)-1)  )  =  YES; 

DISPLAY  0KLVL2; 

OPTIONS  DECIMALS  =8; 

PARAMETER 

P(I,J)  negative  binomial  paramtr 

Q(I,J)  one  minus  the  p  paramtr 

PP(I,J)  p  paramtr  raised  to  power  of  nobs  plus  alpha 

TWT(I,J.K)  total  snipping  wt  of  k  spares  of  component  j 

C(I,J,K)  cumulative  density  for  negative  binomial 

FCI,J,K)  natural  logarithm  of  cumulative  density 

PTMPC I , J , K, L )   temporary  paramtr  needed  to  sum  cdf  recursively; 

P(I, J)$(OK(I, J)  EQ  1)  = 

(TXPOSECI)  +  BETACJ))/((EQDENSCI)5(DURMIS(I))  +  TXPOSE(I)  +  BETA(J)); 

Q(I,J)$(OK(I,J)  EQ  1)  =  1  -  P(I,J); 

PP(I,J)$(OK(I, J)  EQ  1)  =  P(I,J)  XX  (NOBS(J)  +  ALPHA(J)); 

TWT(I,J,K)$0KLVL1(I,J,K)  =  SHIPWT(J)  x  (ORD(K)  -  1); 

PTMP(I,J,K,"DO")$OKLVLl(I,J,k)  =  1; 

X   The  following  LOOP  function  is  used  to  calculate  the  cdf  from  the 
X   negative  binomial  pdf  recursively. 
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LOOPCL,  PTMPCI, J,K,L+1)$0KLVL2(I, J,K,L)  = 

QCI, J)X((N0BS(J)+ALPHA(J)+0RDCL)-1)/C0RD(L)))*PTMP(I,J,K,L)  ); 

C(I,J,K)$0KLVL1(I,J.K)  = 
SUMCL$((0RD(L)-1)  LE  ( LMT( I , J )+ORD( K)-l) ) ,  PTMP( I , J , K, L) )XPP( I , J ) ; 

F(I, J,K)$0KLVL1(I, J,K)  =  LOG(  C(I,J,K)  ); 

DISPLAY  F; 

VARIABLES 

X(I,J,K)   1  if  stock  k  spares  for  component  j  of  end  item  i 
MOE        natural  log  of  probability  of  mission  completion; 

BINARY  VARIABLE  X; 

EQUATIONS 

WTCONST  weight  constraint  equation 

SELECT(I,J)  select  exactly  1  of  the  possible  stockage  levels 

MOEDEF  measure  of  effectiveness  definition; 

WTCONST..   SUM((I,J,K)$0KLVL1(I,J,K),  TWT( I , J , K)XX( I , J , K) )  =L=  WTCAP; 

SELECTCI, J)$C  OK(I,J)  EQ  1  ).. 

SUMC  K$0KLVL1(I,J,K),  X(I,J,K)  )  =E=  1; 

MOEDEF..   SUM((I,J,K)$0KLVL1(I,J,K),  F(  I ,  J  ,  K)5«X(  I ,  J  ,  K)  )  =E=MOE; 

MODEL  STOCK  /ALL/; 

OPTIONS  OPTCR  =  0.0015,  OPTCA  =  0; 

SOLVE  STOCK  USING  MIP  MAXIMIZING  MOE; 

PSUC  =  EXP(  MOE.L  ) 

DISPLAY  X.L, MOE.L, PSUC; 
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